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

A panomics-driven framework for the improvement of major food legume crops: advances, challenges, and future prospects

Food legume crops, including common bean, faba bean, mungbean, cowpea, chickpea, and pea, have long served as vital sources of energy, protein, and minerals worldwide, both as grains and vegetables. Advancements in high-throughput phenotyping, next-generation sequencing, transcriptomics, proteomics, and metabolomics have significantly expanded genomic resources for food legumes, ushering research into the panomics era. Despite their nutritional and agronomic importance, food legumes still face constraints in yield potential and genetic improvement due to limited genomic resources, complex inheritance patterns, and insufficient exploration of key traits, such as quality and stress resistance. This highlights the need for continued efforts to comprehensively dissect the phenome, genome, and regulome of these crops. This review summarizes recent advances in technological innovations and multi-omics applications in food legumes research and improvement. Given the critical role of germplasm resources and the challenges in applying phenomics to food legumes-such as complex trait architecture and limited standardized methodologies-we first address these foundational areas. We then discuss recent gene discoveries associated with yield stability, seed composition, and stress tolerance and their potential as breeding targets. Considering the growing role of genetic engineering, we provide an update on gene-editing applications in legumes, particularly CRISPR-based approaches for trait enhancement. We advocate for integrating chemical and biochemical signatures of cells ('molecular phenomics') with genetic mapping to accelerate gene discovery. We anticipate that combining panomics approaches with advanced breeding technologies will accelerate genetic gains in food legumes, enhancing their productivity, resilience, and contribution to sustainable global food security.

Glycosyltransferases: Pioneering roles in agriculture and medicine

Glycosyltransferases (GTs) belong to a diverse family of enzymes that catalyze the transfer of sugar moieties from activated donor sugars to specific acceptors, thus playing a crucial roles in various biological processes. This review explores the pioneering roles of uridine diphosphate-dependent GTs (UGTs), which use uridine diphosphate glucose as donors. UGTs have also been extensively studied in agricultural and medical fields, emphasizing their potential to revolutionize these sectors. In the agricultural sector, the genetic engineering of UGTs has demonstrate potential in developing crops with enhanced stress tolerance, regulated plant development, and increased resistance to pests and diseases. These advancements not only contribute to sustainable farming practices but also address global food security challenges by facilitating the production of more resilient plant varieties. Furthermore, UGTs facilitate the synthesis of complex carbohydrates and glycoconjugates in plants, which are critical for developing drugs and therapeutic strategies targeting various ailments, including cancer and infectious diseases. Thus, this review explored the functions and synthesis methods of flavonoid glycosides, terpenoid glycosides, and polyketosides in detail. Moreover, owing to the functional diversity of UGTs, numerous research methods were reviewed, and novel, more valuable UGTs will be obtained. In summary, this study synthesizes the current research findings and discusses future perspectives to underscore the transgenic technology and synthetic biological impact of UGTs on agriculture and medicine and bridge the gap between fundamental science and practical applications.

An RNase H-Like gene complements resistance to Bean common mosaic necrosis virus in Phaseolus vulgaris

Bean common mosaic virus and Bean common mosaic necrosis virus (BCMNV) are related positive-sense RNA potyviruses that limit the production of common bean (Phaseolus vulgaris L.) worldwide. Potyviruses infect other legume species, such as Glycine max, which can serve as a source of orthologous resistance genes. The bc-1 on chromosome Pv03 in P. vulgaris and Rsv4 on Gm02 in G. max are syntenic gene regions that confer resistance to Potyviruses. Receptor-like kinases underlie both gene regions, and a linked RNase H-Like protein, which enhances potyvirus resistance, was recently associated with Rsv4. We sought to determine if RNase H-Like encoding genes are similarly located in the bc-1 region and enhance resistance to BCMNV. Synteny analysis between P. vulgaris and G. max revealed RNase H-Like genes near the bc-1 locus. Positional cloning among reference genomes and select genotypes, along with linkage mapping in recombinant inbred line and F2 populations, identified the RNase H-Like 1 gene, which enhanced resistance to BCMNV NL-3 strain when combined with I, bc-1, and bc-u or with bc-1 and bc-u genes by reducing systemic spread of susceptible symptoms in non-inoculated leaves. A single nucleotide polymorphism marker, G03_4166082, was developed to track the resistant and susceptible alleles for RNase H-Like 1 in breeding programs. Overall, this study advances the understanding of the complex mechanisms underlying BCMNV resistance in common bean.

Evolution of NAC transcription factors from early land plants to domesticated crops

NAC [NO APICAL MERISTEM (NAM), ARABIDOPSIS TRANSCRIPTION ACTIVATOR FACTOR 1/2 (ATAF1/2), and CUP-SHAPED COTYLEDON (CUC2)] transcription factors are key regulators of plant growth, development, and stress responses but were also crucial players during land plant adaptation and crop domestication. Using representative members of green algae, bryophytes, lycophytes, gymnosperms, and angiosperms, we expanded the evolutionary history of NAC transcription factors to unveil the relationships among members of this gene family. We found a massive increase in the number of NAC transcription factors from green algae to lycophytes and an even larger increase in flowering plants. Many of the NAC clades arose later during evolution since we found eudicot- and monocot-specific clades. Cis-elements analysis in NAC promoters showed the presence of abiotic and biotic stress as well as hormonal response elements, which indicate the ancestral function of NAC transcription factor genes in response to environmental stimuli and in plant development. At the transcriptional level, the expression of NAC transcription factors was low or absent in male reproduction, particularly mature pollen, across the plant kingdom. We also identified NAC genes with conserved expression patterns in response to heat stress in Marchantia polymorpha and Oryza sativa. Our study provides further evidence that transcriptional mechanisms associated with stress responses and development emerged early during plant land adaptation and are still conserved in flowering plants and domesticated crops.

Common bean pan-genome reveals abundant variation patterns and relationships of stress response genes and pathways

Long-term geographical isolation and the different directions of domestication can cause a large number of genome variations. Population genetic analysis based on a single reference genome cannot capture all the variation information. Pan-genome construction is an effective way to overcome this problem. Resequencing data from 683 common bean landraces and breeding lines provided a pan-genome construction data resource. For the first time, for common bean pan-genome construction, 305 Mb non-reference contigs and 10,452 novel genes were identified. Among these new genes, 373 resistance gene analogs containing 372 variable genes were identified and used to narrow down the candidate genes in Pseudomonas syringae pv. phaseolicola resistance quantitative trait locus interval of the common bean. Transcriptome analysis of multiple biotic and abiotic stresses reveals that gene expression patterns are organ-, stress-, and gene conservation-specific. Core and shell genes may be co-expressed in all samples and may have functional complementarity to maintain the stability of plant growth. Within pathways, 8990 and 30,272 mutual exclusivity and co-occurrence gene presence-absence variations (PAVs) were discovered respectively, providing further insights into the functional complementarity of genes. In conclusion, our study provides a comprehensive genome resource, which will be useful for further common bean breeding and study.

Genomics reveal population structure, genetic diversity and evolutionary history of Phyllostachys edulis (moso bamboo) in global natural distribution

Introduction

Moso bamboo (Phyllostachys edulis) is widespread in natural forests over large areas in China.

Methods

Here we collected 193 individuals of moso bamboo from 37 natural populations in China's distribution area. Genotyping by sequencing (GBS) was employed to elucidate the genetic diversity, genetic structure, selection pressure, history and adaptive distribution prediction of moso bamboo.

Results

The results revealed that the moso bamboo in China can be divided into central α, eastern β and southern γ subpopulations, with the α-subpopulation presumed to be the origin center. Notably, the genetic diversity of moso bamboo populations were relatively low, and the heterozygotes were excess. At the subpopulation level, the genetic diversity of α-subpopulation was the highest and that of β-subpopulation was the lowest. Analysis of population selection pressure during the transmission of moso bamboo indicated significant genetic differences among subpopulations. Furthermore, 3681 genes related to adaptability, stress resistance, photosynthesis, and hormone were identified from the selected regions. Four SNP markers developed and validated. Based on the population dynamics history and distribution simulation, we found that the distribution of moso bamboo has been influenced by the climate change in geological history.

Discussion

These findings hold significant implications for enhancing our genetic comprehension of bamboo populations and exploring germplasm resources.

Evaluation of the Effect of Light Color on Albumins and Globulins Content During Bean Germination

The effect of different light colors on the concentration of albumins and globulins during black bean (Phaseolus vulgaris L.) germination was evaluated with an RGB LED lighting system. This study aimed to determine how light of different spectral composition influences protein content across different germination stages. Black bean seeds were germinated under six different LED light sources (red, green, blue, white, violet, and cyan), and protein fractions were quantified by Bradford's method. The results showed that blue and cyan light increased the concentration of albumins. Blue, white, and cyan light for globulins increased the concentration compared to germination under solar conditions for both fractions. The positive correlation between these protein fractions under specific wavelengths suggests a metabolic adaptation to light exposure. These findings highlight the potential of controlled lighting conditions to enhance the nutritional quality of germinated beans, supporting their application as functional food ingredients. Additionally, this study underscores the importance of photobiological modulation in seed germination, opening new possibilities for optimizing plant-based protein sources. Future research could explore the mechanisms behind these protein variations and their implications for food production and nutrition.

The unexplored diversity of rough-seeded lupins provides rich genomic resources and insights into lupin evolution

Lupin crops provide nutritious seeds as an excellent source of dietary protein. However, extensive genomic resources are needed for crop improvement, focusing on key traits such as nutritional value and climate resiliency, to ensure global food security based on sustainable and healthy diets for all. Such resources can be derived either from related lupin species or crop wild relatives, which represent a large and untapped source of genetic variation for crop improvement. Here, we report genome assemblies of the cross-compatible species Lupinus cosentinii (Mediterranean) and its pan-Saharan wild relative L. digitatus, which are well adapted to drought-prone environments and partially domesticated. We show that both species are tetraploids, and their repetitive DNA content differs considerably from that of the main lupin crops L. angustifolius and L. albus. We present the complex evolutionary process within the rough-seeded lupins as a species-based model involving polyploidization and rediploidization. Our data also provide the foundation for a systematic analysis of genomic diversity among lupin species to promote their exploitation for crop improvement and sustainable agriculture.

Expression profiling and transcriptional regulation of the SRS transcription factor gene family of common bean (Phaseolus vulgaris) in symbiosis with Rhizobium etli

The SRS/STY transcription factors from the model legumes: Lotus japonicus and Medicago truncatula, are part of regulatory networks that play relevant roles for nodule development during the N-fixing symbiosis with rhizobia. In this work we analyzed the participation of the PvSRS transcription factors from common bean (Phaseolus vulgaris), a most important legume crop, in the symbiosis with Rhizobium etli. Our phylogenetic analysis of SRS TFs across five plant species, including four legumes and Arabidopsis thaliana, identified clades that group SRS proteins that are highly expressed in legume nodules and in Arabidopsis roots. A qRT-PCR expression analysis of the 10 PvSRS in root/nodule of inoculated plants, revealed that all the PvSRS genes are expressed at different stages of the symbiosis, albeit at different levels. Based on what is known for L. japonicus, we demonstrated that the PvSRS10 gene -with highest expression during symbiosis- is transcriptionally activated by NF-Y transcription factor, thus indicating its participation in the NIN-NF-Y regulatory cascade. Based on our previous work about the relevant role of members from the MADS-domain/AGL transcription factors as regulators of the N-fixing symbiosis, in this work we demonstrated the transcriptional regulation of PvSRS10 by the MADS-TF PvFUL-like. Analysis of protein-protein interaction networks predicted thatPvSRS5 and PvSRS6 interact with proteins involved in transcriptional regulation and the auxin-activated signaling pathway. The regulatory mechanisms of PvSRS TF in common bean symbiosis may be related to auxin biosynthesis regulation, that is essential for determinate nodules development. Our study highlights the role of PvSRS TF in the N-/fixing symbiosis, a relevant process for sustainable agriculture.

Whole plant transpiration responses of common bean (Phaseolus vulgaris L.) to drying soil: Water channels and transcription factors

Common bean (Phaseolus vulgaris L.) is the main legume crop for direct human consumption worldwide. Among abiotic factors affecting common bean, drought is the most limiting. This study aimed at characterizing genetic variability and architecture of transpiration, stomatal regulation and whole plant water use within the Mesoamerican germplasm. A critical fraction of transpirable soil water (FTSWc) was estimated as the inflection point at which NTR starts decreasing linearly. Genome-wide association (GWA) analyses for mean NTR and FTSWc were performed. High variation on mean NTR and FTSWc was found among genotypes. Unreported genomic signals controlling the variation of these traits were identified on Pv01 and Pv07 some located in intergenic, intronic and exonic regions. A set of novel candidate genes and putative regulatory elements located in these QTL were identified. Some of the genes have been previously reported to be involved in abiotic tolerance in model species, including some of the five transcription factors (TF) identified. Four candidate genes, one with potential water transportation activity and three TFs were validated. The gene Phvul.001G108800, an aquaporin SIP2-1 related gene, showed water channel activity through oocyte water assays. Mutant Arabidopsis thaliana (Ath) lines for the homologous genes of common bean were evaluated in transpiration experiments. Two of the three evaluated TFs, UPBEAT1 and C2H2-type ZN finger protein, were involved in the control of transpiration responses to drying soil. Our results provide evidence of novel genes to accelerate the drought tolerance improvement in the crop and study the physiological basis of drought response in plants.

Genotype-Specific Responses of Common Bean to Meloidogyne incognita

The root-knot nematode Meloidogyne incognita causes large galls on roots, interfering with the flow of water and nutrients to the plant. In the common bean, no major resistance (R) genes have been described. Instead, resistance is controlled by multiple genes, which have not proved effective so far. An RNA sequencing approach was used to identify genes involved in common bean response to M. incognita at the stages of nematode invasion and root galling. When comparing infected and uninfected treatments of a moderately resistant (MR) and a susceptible (S) genotype, several genes were identified as differentially expressed. Their functional annotation indicated that both genotypes underwent complex transcriptional reprogramming from early to later periods of the interaction, but defense-related genes were mostly upregulated in the MR genotype. At the early stage, a large set of genes was activated in both genotypes, including those involved in cell wall organization, signaling, hormonal pathways, transcription factors, oxidative stress, and putative resistance gene analogs. Later, most of the previously activated defense mechanisms were no longer expressed in the S genotype. There was an increased expression of genes encoding proteins involved in hormonal signaling pathways (salicylic acid and gibberellin-related), protein kinases, transcription factors, and oxidative stress in the MR genotype. However, a decreased expression of genes involved in signaling mediated by calcium and oxidative stress occurred in the S genotype, indicating susceptibility. The repertoire of genes identified herein will facilitate research in plant-nematode interactions, with possible applications for the improvement of the common bean.

Crop assemblages in the highlands of Odisha: toward an ethnoarchaeobotanical perspective of crop choices

The highlands of Odisha are home to numerous indigenous communities, known as Adivasi, who traditionally cultivate various species of millets and rice in systems of shifting cultivation and permanent upland cultivation. National and international agricultural schemes have had notable impact on the crops cultivated in the highlands, increasing the consumption of (modern varieties of) rice at the expense of millets. Recent initiatives, such as the Odisha Millet Mission, are attempting to reverse the dependence on modern rice and reintroduce millets, in the highlands as well as in the lowland urban areas. The archaeobotanical records of South Asia and beyond illustrate that crop assemblages have always been dynamic, and understanding the choices that underlie these dynamics can contribute insights into the emergence of contemporary crop assemblages and related foodways, as well as their future trajectories. This contribution seeks to explore how crop choices are reflected in current crop assemblages in the highlands of Odisha and to examine how anthropological insights can contribute to an understanding of crop choices in the past and vice versa.This article is part of the theme issue 'Unravelling domestication: multi-disciplinary perspectives on human and non-human relationships in the past, present and future'.

The genomic landscape of gene-level structural variations in Japanese and global soybean Glycine max cultivars

Japanese soybeans are traditionally bred to produce soy foods such as tofu, miso and boiled soybeans. Here, to investigate their distinctive genomic features, including genomic structural variations (SVs), we constructed 11 nanopore-based genome references for Japanese and other soybean lines. Our assembly-based comparative method, designated 'Asm2sv', identified gene-level SVs comprehensively, enabling pangenome analysis of 462 worldwide cultivars and varieties. Based on these, we identified selective sweeps between Japanese and US soybeans, one of which was the pod-shattering resistance gene PDH1. Genome-wide association studies further identified several quantitative trait loci that accounted for large-seed phenotypes of Japanese soybean lines, some of which were also close to regions of the selective sweeps, including PDH1. Notably, specific combinations of alleles, including SVs, were found to increase the seed size of some Japanese landraces. In addition to the differences in cultivation environments, distinct food processing usages might result in changes in Japanese soybean genomes.

Genome-wide Identification, Characterization, and Expression Analysis of NHX Genes in Phaseolus vulgaris L. under Salt Stress: An In Silico Approach

BACKGROUND

Climate change is among the major triggering agents of abiotic stresses (e.g., saline stress), culminating in a vulnerability of common bean production systems. In recent decades, important research has identified and characterized genes that can mitigate the adverse effects caused by salt stress; among them, the Na+/H+ antiporters (NHXs) gene stands out. The NHX genes are widely distributed in all organisms and play significant roles in osmotic regulation in plants under salt stress conditions. Genome-wide identification of NHX genes has been performed in several plant species but not in Phaseolus vulgaris L.

METHODS

This study aimed to identify and characterize NHX genes in P. vulgaris L. using a genome-wide analysis approach conducted in silico. The common bean genome revealed nine putative PvNHX genes, and their subcellular localization, phylogenetic relationship, cis-regulatory elements, conserved motifs identification, chromosomal location, expression patterns, and interaction networks were analyzed.

RESULTS

Promoter analysis suggested that PvNHX genes shared hormone-related elements and were light-responsive and stress-responsive. Seven PvNHX genes were under the regulation of five microRNA (miRNA) families. RNA-seq analysis revealed that most PvNHX genes were expressed in response to salt stress. Currently, the most assertive strategy to confront these adversities is to use the information generated by sequencing plants to identify candidate genes that can be introgressed to improve programs in producing resilient cultures.

CONCLUSION

These results can provide valuable information for future studies on the functional mechanism of PvNHX genes in common beans in response to salt stress.

GWAS-assisted and multitrait genomic prediction for improvement of seed yield and canning quality traits in a black bean breeding panel

In recent years, black beans (Phaseolus vulgaris L.) have gained popularity in the United States, with improved seed yield and canning quality being critical traits for new cultivars. Achieving genetic gains in these traits is often challenging due to negative trait associations and the need for specialized equipment and trained sensory panels for evaluation. This study investigates the integration of genomics and phenomics to enhance selection accuracy for these complex traits. We evaluated the prediction accuracy of single-trait (ST) and multitrait (MT) genomic prediction (GP) models, incorporating near-infrared spectroscopy (NIRS) data and markers identified through genome-wide association studies (GWAS). The models demonstrated moderate prediction accuracies for yield and canning appearance (App) and high accuracies for color retention. No significant differences were found between ST and MT models within the same breeding cycle. However, across breeding cycles, MT models outperformed ST models by up to 45 and 63% for canning App and seed yield, respectively. Interestingly, incorporating significant SNP markers identified by GWAS and NIRS data into the models tended to decrease prediction accuracy both within and between breeding cycles. As genotypes from the new breeding cycle were included, the models' prediction accuracy generally increased. Our findings underscore the potential of MT models to enhance the prediction of complex traits such as seed yield and canning quality in dry beans and highlight the importance of continually updating the training dataset for effective GP implementation in dry bean breeding.

A special short-wing petal faba genome and genetic dissection of floral and yield-related traits accelerate breeding and improvement of faba bean

BACKGROUND

A comprehensive study of the genome and genetics of superior germplasms is fundamental for crop improvement. As a widely adapted protein crop with high yield potential, the improvement in breeding and development of the seeds industry of faba bean have been greatly hindered by its giant genome size and high outcrossing rate.

RESULTS

To fully explore the genomic diversity and genetic basis of important agronomic traits, we first generate a de novo genome assembly and perform annotation of a special short-wing petal faba bean germplasm (VF8137) exhibiting a low outcrossing rate. Comparative genome and pan-genome analyses reveal the genome evolution characteristics and unique pan-genes among the three different faba bean genomes. In addition, the genome diversity of 558 accessions of faba bean germplasm reveals three distinct genetic groups and remarkable genetic differences between the southern and northern germplasms. Genome-wide association analysis identifies several candidate genes associated with adaptation- and yield-related traits. We also identify one candidate gene related to short-wing petals by combining quantitative trait locus mapping and bulked segregant analysis. We further elucidate its function through multiple lines of evidence from functional annotation, sequence variation, expression differences, and protein structure variation.

CONCLUSIONS

Our study provides new insights into the genome evolution of Leguminosae and the genomic diversity of faba bean. It offers valuable genomic and genetic resources for breeding and improvement of faba bean.

A multidisciplinary and integrative review of the structural genome and epigenome of Capsicum L. species

MAIN CONCLUSION

We revised and integrated the genomic and epigenomic data into a comparative Capsicum ideogram, evidencing the advances and future perspectives. Capsicum L. (Solanaceae) genome has been characterized concerning karyotype, nuclear and chromosomal genome size, genome sequencing and physical mapping. In addition, the epigenome has been investigated, showing chromosomal distribution of epimarks in histone amino acids. Genetic and epigenetic discoveries have given light to understanding the structure and organization of the Capsicum "omics". In addition, interspecific and intraspecific similarities and diversities have been identified, characterized and compared in taxonomic and evolutive scenarios. The journey through Capsicum studies allows us to know the 2n = 2x = 24 and 2n = 2x = 26 chromosome numbers, as well as the relatively homomorphic karyotype, and the 1C chromosomal DNA content. In addition, Capsicum "omics" diversity has mainly been evidenced from the nuclear 1C value, as well as from repeatome composition and mapping. Like this, Capsicum provides several opportunities for "omics", ecological, agronomic and conservation approaches, as well as subjects that can be used at different levels of education. In this context, we revisit and integrate Capsicum data about the genome size, karyotype, sequencing and cytogenomics, pointing out the progress and impact of this knowledge in taxonomic, evolutive and agronomic contexts. We also noticed gaps, which can be a focus of further studies. From this multidisciplinary and integrative review, we intend to show the beauty and intrigue of the Capsicum genome and epigenome, as well as the outcomes of these similarities and differences.

Fine mapping of the unique Ur-11 gene conferring broad resistance to the rust pathogen of common bean

KEY MESSAGE

Fine mapping positioned the Ur-11 rust resistance gene in common bean to a narrow 9 kb genomic region and enabled the development of a KASP marker tightly linked to Ur-11 for use in gene pyramiding to achieve durable rust resistance. The extensive virulence diversity of the fungal pathogen Uromyces appendiculatus threatens common bean (Phaseolus vulgaris) production. The Ur-11 gene present in the Guatemalan common bean accession PI 181996 conferred resistance to 89 of 90 virulent races of U. appendiculatus. We describe here the fine mapping of Ur-11 and the development and validation of a DNA marker tightly linked to Ur-11. An F2 population from the cross between the susceptible Pinto 114 with the resistant PI 181996 was inoculated with four races of U. appendiculatus. This study established that the rust resistance in PI 181996 was conferred by Ur-11. We then fine mapped Ur-11 using F2 plants and F2:3 families, high-throughput SNP genotyping, SSRs and KASPs marker development, whole-genome sequencing, and local haplotype analysis. Ur-11 was positioned in a narrow 9.01 Kb genomic region on chromosome Pv11 flanked by KASP markers SS322 and SS375. This genomic region included a candidate gene encoding a nucleotide-binding site and leucine rich-repeat domain with pathogen resistance functions. The validation of the SS322 KASP marker was performed on a panel of 206 diverse common bean cultivars that were inoculated with four races of U. appendiculatus. The SS322 marker was 97.5% accurate in identifying the presence of Ur-11 in common bean plants. These results suggest that S322 will be a highly effective molecular marker for the development of common bean cultivars with Ur-11 alone and combining Ur-11 with other rust resistance genes that would confer broad and durable resistance to the hypervirulent bean rust pathogen.

The Genome of the Lima Bean Variety Baiyu Bean Highlights Its Evolutionary Characteristics

The baiyu bean (Phaseolus lunatus), also known as the lima bean, is a plant belonging to the Fabaceae family, has a long and distinguished history of cultivation in China and is a highly regarded local variety of lima bean. In the current study, we present the reference genome of the baiyu bean variety, which has a scaffold N50 length of 47.545 Mb. A comparative genomic analysis was conducted using genomes of seven legume species, and the results demonstrated that 1564 and 1275 genes of baiyu bean exhibited expansion and contraction, respectively. Moreover, 543 genes were identified as exclusive to the baiyu bean. The analysis of adaptive evolution genes revealed the presence of 61 genes under adaptive evolution between P. lunatus and the common bean P. vulgaris. An examination of the branch model revealed the presence of five genes undergoing adaptive evolution in the P. lunatus branch. Additionally, the evolutionary selective pressure acting on other branches of legume plants was analyzed. A comprehensive analysis of structural variations (SVs) between the baiyu bean and G27455 genome was conducted, resulting in the identification of 5549 SVs. Among these, 333 genes were identified as high-impact SV genes. The acquisition of the genome sequence of this excellent variety will facilitate the exploration and utilization of its characteristics, providing a foundation for the genetic improvement of the lima bean.

Predicting the role of β-GAL genes in bean under abiotic stress and genome-wide characterization of β-GAL gene family members

Β-Gals are a subgroup of the glycoside hydrolase (GH) family of enzymes, which possess the Glyco_hydro_35 (GH35) domain. Although studies have been conducted on the β-Gal gene family in numerous plant species, no such research has been conducted on beans. The purpose of this study was to determine the gene expression levels of β-Gal genes in the leaf tissue of P. vulgaris under salt and drought stress using quantitative real-time polymerase chain reaction (qRT-PCR) and to perform a comprehensive analysis of β-Gal gene family members using bioinformatics tools. In the bean genome, 25 Pvul-βGAL proteins with amino acid numbers ranging from 291 to 1119, molecular weights from 32.94 to 126.56 kDa, and isoelectric points from 5.46 to 9.08 were identified. Both segmental and tandem duplication have occurred in β-Gal genes in the bean genome, and Pvul-BGAL genes have been subject to negative selection in the evolutionary process. For a deeper comprehension of the evolutionary proximity of Pvul-BGAL genes, a phylogenetic tree and synteny map were drawn together with Arabidopsis thaliana and Glycine max β-Gal genes. The expression profiles of β-Gal genes in different tissues of the bean were determined in silico. In addition, the expression profiles of β-Gal genes in the leaves of bean plants subjected to drought and salt stress were analyzed, and the role of β-Gal genes in salt and drought stress was estimated. In this study, the role of β-Gal gene family in abiotic stress response and the characterization of β-Gal genes in beans were determined for the first time and will provide a basis for future functional genomics studies.

Chromosome-level de novo genome unveils the evolution of Gleditsia sinensis and thorns development

Gleditsia sinensis Lam. (G. sinensis) as an important species within the Leguminosae family, has been utilized in Chinese medicine for centuries, and its thorns serve as a chief medicinal ingredient. The absence of a comprehensive genome database has hindered its in-depth research. In this investigation, a chromosome-level de novo genome assembly of G. sinensis 'Yulin No.1' was achieved, which harbors a 786.13 Mb sized genome with 36,408 protein-coding genes and experiences two WGD events. The comparative and evolutionary analysis unveiled the close phylogenetic relationship between G. sinensis and eight other Leguminosae species. The WGCNA and gene family analysis further indicated that GsinMYB was involved in the development of thorns. This investigation offered a high-level genome of G. sinensis, facilitating comparisons in Leguminosae species evolution and functional elucidation. It also provided key insights for further research on the molecular regulation mechanisms of thorn development in plants and the molecular breeding of G. sinensis.

Mapping resistance to Sclerotinia white mold in two pinto bean recombinant inbred line populations

White mold, caused by the fungus Sclerotinia sclerotiorum (Lib.) de Bary, is a devastating disease affecting common bean (Phaseolus vulgaris L.) production worldwide. Breeding for resistance to white mold is challenging due to its quantitative inheritance and intricate genetic mechanisms. This research aimed to validate and characterize physiological resistance in the pinto dry bean market class through greenhouse straw tests under controlled conditions and field assessments under natural environments. Classical quantitative trait locus (QTL) mapping and Khufu de novo QTL-seq were employed to detect and narrow QTL intervals and identify candidate genes associated with white mold resistance in two pinto bean recombinant inbred line populations, PT9-5-6/USPT-WM-12 (P2) and PT12-37/VCP-13 (P3). Eleven QTL, five in P2 and six in P3, conditioning white mold resistance were identified. New QTL were discovered including WM1.4 and WM11.5 in P2, and WM1.5 and WM7.7 in P3. Existing major-effect QTL were validated: WM5.4 (34%-phenotypic variation explained) and WM7.4 (20%) in straw tests, and WM2.2 (15%) and WM3.1 (27%) under field conditions. QTL for avoidance traits such as resistance to lodging and late maturity overlapped WM2.2 in P2 and WM1.5, WM3.1, WM5.4, and WM7.7 in P3. WM5.4 (Pv05: 7.0-38.7 Mb) was associated with a large Phaseolus coccineus L. genome introgression in the resistant parent VCP-13. These findings offer narrowed genomic intervals and putative candidate genes for marker-assisted selection targeting white mold resistance improvement in pinto beans.

Profiling conserved transcription factor binding motifs in Phaseolus vulgaris through comparative genomics

Common bean (Phaseolus vulgaris), a staple food in Latin America and Africa, serves as a vital source of energy, protein, and essential minerals for millions of people. However, genomics knowledge that breeders could leverage for improvement of this crop is scarce. We have developed and validated a comparative genomics approach to predict conserved transcription factor binding sites (TFBS) in common bean and studied gene regulatory networks. We analyzed promoter regions and identified TFBS for 12,631 bean genes with an average of 6 conserved motifs per gene. Moreover, we discovered a statistically significant relationship between the number of conserved motifs and amount of available experimental evidence of gene regulation. Notably, ERF, MYB, and bHLH transcription factor families dominated conserved motifs, with implications for starch biosynthesis regulation. Furthermore, we provide gene regulatory data as a resource that can be interrogated for the regulatory landscape of any set of genes. Our results underscore the significance of TFBS conservation in legumes and aligns with the notion that core genes often exhibit a more conserved regulatory makeup. The study demonstrates the effectiveness of a comparative genomics approach for addressing genome information gaps in non-model organisms and provides valuable insights into the regulatory networks governing starch biosynthesis genes that can support crop improvement programs.

The evolutionary history of the common bean (Phaseolus vulgaris) revealed by chloroplast and nuclear genomes analysis

KEY MESSAGE

The origin of common bean was investigated throughout chloroplast and nuclear WGS data considering recombination events. Our results support the Mesoamerican origin of common bean. The remarkable evolutionary history of the common bean (Phaseolus vulgaris L.) has led to the emergence of three wild main gene pools corresponding to three different eco-geographical areas: Mesoamerica, the Andes and northern Peru/Ecuador. Recent works proposed novel scenarios, and the northern Peru/Ecuador population has been described as a new species called P. debouckii, rekindling the debate about the origin of P. vulgaris. Here we shed light on the origin of P. vulgaris by analyzing the chloroplast and nuclear genomes of a large varietal collection representing the entire geographical distribution of wild forms including a large collection of Mesoamerican and Andean individuals. We assembled 37 chloroplast genomes de novo and used them to construct a time frame for the divergence of the genotypes under investigation, revealing that the separation of the Mesoamerican and northern Peru/Ecuador gene pools occurred ~ 0.15 Mya. Our results clearly support a Mesoamerican origin of the common bean and reject the recent P. deboukii hypothesis. These results also imply two independent migratory events from Mesoamerica to the North and South Andes, probably facilitated by birds. Our work represents a paradigmatic example of the importance of taking into account the genetic rearrangements produced by recombination when investigating phylogeny and of the analysis of wild forms when studying the evolutionary history of a crop species.

Genomic insights into genetic diversity and seed coat color change in common bean composite populations

Introduction

The color of the seed coat of common bean (Phaseolus vulgaris L.) is an important trait influencing marketability and consumer preferences. An understanding of the genetic mechanisms underlying seed coat color variation can aid in breeding programs aimed at improving esthetic and agronomic traits. This study investigates the genetic diversity and molecular mechanisms associated with seed coat color change in composite bean populations through phenotypic analysis and whole genome sequencing (WGS).

Methods

Four composite populations and two standard varieties of common bean were cultivated over a two-year period and seed coat color and morphological traits were assessed. WGS was performed on 19 phenotypes and yielded 427 GB of data with an average sequencing depth of 30×. More than 8.6 million high-confidence single nucleotide polymorphisms (SNPs) were identified. Genetic diversity metrics such as nucleotide diversity (π), observed heterozygosity (Ho), expected heterozygosity (He) and allelic richness (Ar) were calculated. Population structure was analyzed using Fst, principal component analysis (PCA) and clustering. Cross-population statistics (XP-CLR and XP-EHH) were used to identify selection signals associated with seed coat color change. Gene Ontology (GO) and KEGG enrichment analyzes were performed for candidate genomic regions.

Results

Phenotypic analysis revealed significant differences in seed coat color among the four composite populations, with notable changes among years. The populations exhibited different growth habits and plant types, especially KIS_Amand and SRGB_00366, which showed the highest phenotypic diversity in seed coat color. WGS identified 8.6 million SNPs, with chromosomes 4 and 1 having the highest SNP density (11% each), while chromosomes 3 and 6 had the lowest. KIS_Amand had the highest genetic diversity (π = 0.222, Ar = 1.380) and SRGB_00189 the lowest (π = 0.067, Ar = 1.327). SRGB_00366 showed moderate genetic diversity (π = 0.173, Ar = 1.338) and INCBN_03048 showed medium diversity (π = 0.124, Ar = 1.047). The Fst values indicated a strong genetic differentiation, especially between the two standard varieties ETNA and Golden_Gate (Fst = 0.704) and the composite populations. Selective sweep analysis with XP-CLR and XP-EHH identified 118 significant regions associated with seed coat color change, with most regions located on chromosomes 4, 9, 10 and 11. Phosphatidylinositol signaling pathways were highly enriched in candidate regions, indicating that cellular transport mechanisms play a critical role in seed coat pigmentation. Key GO terms included phosphatidylinositol-biphosphate binding, exocytosis, and vesicle-mediated transport, suggesting a link between cellular transport and pigment deposition in the seed coat.

Discussion

The study demonstrates significant genetic diversity within and among common bean composite populations, with KIS_Amand and SRGB_00366 exhibiting the highest phenotypic and genetic variability. The identification of selective sweeps and the enrichment of phosphatidylinositol-related pathways provide new insights into the molecular mechanisms controlling seed coat color variation. The strong genetic differentiation between standard varieties and composite populations highlights the role of selective breeding in shaping the genetic landscape of common bean. The results suggest that variation in seed coat color is controlled by both regulatory and structural genetic changes, providing valuable information for breeding programs.

Conclusion

This study provides a detailed analysis of the genetic architecture of seed coat color variation in common bean. The identification of key genomic regions and pathways associated with seed pigmentation improves our understanding of the complex genetic interactions underlying this trait. These results provide valuable genomic resources for future breeding efforts aimed at improving seed color and other important traits in common bean.

High-throughput sequencing revealed the symptomatic common bean (Phaseolus vulgaris L.) virome in Iran

Common bean (Phaseolus vulgaris L.) is a crop rich in protein, minerals, and starch. Viruses are a significant limiting factor in increasing the production of legumes, particularly common beans. Accurate and timely detection of plant viruses is essential for minimizing disease damage and ensuring food security. To investigate common bean field viruses in Iran, 300 samples of common bean plants showing viral-like symptoms were collected over 3 years, 2020, 2021, and 2022. This study is the first to use total RNA-seq for a virome analysis of common beans in Iran. The results of the total RNA-seq indicated that the common bean samples were infected with sesame curly top virus (SeCTV), beet curly top Iran virus (BCTIV), tomato leaf curl Palampur virus (ToLCPalV), cucumber mosaic virus (CMV), bean common mosaic virus (BCMV), phaseolus vulgaris endornavirus 1 (PvEV1) and phaseolus vulgaris endornavirus 2 (PvEV2). This is the first report of PvEV1 and PvEV2 in Iran. Moreover, these findings revealed the presence of SeCTV and ToLCPalV for the first time in the western part of the country. Furthermore, the analysis of the nucleotide sequences and the phylogeny tree obtained from the complete genome of the two BCTIV isolates in this study, compared to other isolates, indicated the presence of a new strain of BCTIV in the common bean fields. During the three-year study, the detection rate of viruses indicated that BCTIV and BCTV were more prevalent in Lorestan province than in Markazi province. The research findings showed that common bean fields in the central and western regions of the country were infected with seven viruses, with DNA viruses being more prevalent in Lorestan province than in Markazi province. This information should be taken into account when developing management strategies and breeding plans.

A dTALE approach demonstrates that induction of common bean OVATE Family Protein 7 promotes resistance to common bacterial blight

Common bacterial blight (CBB) is a devastating seed-transmitted disease of common bean (Phaseolus vulgaris L.), caused by Xanthomonas phaseoli pv. phaseoli and Xanthomonas citri pv. fuscans. The genes responsible for CBB resistance are largely unknown. Moreover, the lack of a reproducible and universal transformation protocol limits the study of genetic traits in common bean. We produced X. phaseoli pv. phaseoli strains expressing artificially designed transcription-activator like effectors (dTALEs) to target 14 candidate genes for resistance to CBB based on previous transcriptomic data. In planta assays in a susceptible common bean genotype showed that induction of PvOFP7, PvAP2-ERF71, or PvExpansinA17 expression by dTALEs resulted in CBB symptom reduction. After PvOFP7 induction, in planta bacterial growth was reduced at early colonization stages, and RNA-seq analysis revealed up-regulation of cell wall formation and primary metabolism, together with major down-regulation of heat shock proteins. Our results demonstrated that PvOFP7 contributes to CBB resistance, and underlined the usefulness of dTALEs for functional validation of genes whose induction impacts Xanthomonas-plant interactions.

CRISPR/Cas9 editing of two adenine phosphoribosyl transferase coding genes reveals the functional specialization of adenine salvage proteins in common bean

Adenine metabolism is important for common bean (Phaseolus vulgaris L.) productivity since this legume uses ureides derived from the oxidation of purine nucleotides as its primary nitrogen storage molecules. Purine nucleotides are produced from de novo synthesis or through salvage pathways. Adenine phosphoribosyl transferase (APRT) is the enzyme dedicated to adenine nucleobase salvage for nucleotide synthesis, but it can also convert active cytokinin bases into their inactive nucleotide forms. In common bean, APRT is encoded by four genes. Gene expression analysis, biochemical properties, and subcellular location indicated functional differences among the common bean APRT isoforms. CRISPR/Cas9 targeted down-regulation of two of the four PvAPRTs followed by metabolomic and physiological analyses of targeted hairy roots revealed that, although the two proteins have redundant functions, PvAPRT1 mostly participated in the salvage of adenine, whereas PvAPRT5 was the predominant form in the regulation of cytokinin homeostasis and stress responses with a high impact in root and nodule growth.

Telomere-to-telomere genome of common bean (Phaseolus vulgaris L., YP4)

BACKGROUND

Common bean is a significant grain legume in human diets. However, the lack of a complete reference genome for common beans has hindered efforts to improve agronomic cultivars.

FINDINGS

Herein, we present the first telomere-to-telomere (T2T) genome assembly of common bean (Phaseolus vulgaris L., YP4) using PacBio High-Fidelity reads, ONT ultra-long sequencing, and Hi-C technologies. The assembly resulted in a genome size of 560.30 Mb with an N50 of 55.11 Mb, exhibiting high completeness and accuracy (BUSCO score: 99.5%, quality value (QV): 54.86). The sequences were anchored into 11 chromosomes, with 20 of 22 telomeres identified, leading to the formation of 9 T2T pseudomolecules. Furthermore, we identified repetitive elements accounting for 61.20% of the genome and predicted 29,925 protein-coding genes. Phylogenetic analysis suggested an estimated divergence time of approximately 11.6 million years ago between P. vulgaris and Vigna angularis. Comparative genome analysis revealed the expanded gene families and variations between YP4 and G19833 associated with defense response.

CONCLUSIONS

The T2T reference genome and genomic insights presented here are crucial for future genetic studies not only in common bean but also in other legumes.

A chromosome-scale genome assembly of mungbean (Vigna radiata)

Background

Mungbean (Vigna radiata) is one of the most socio-economically important leguminous food crops of Asia and a rich source of dietary protein and micronutrients. Understanding its genetic makeup is crucial for genetic improvement and cultivar development.

Methods

In this study, we combined single-tube long-fragment reads (stLFR) sequencing technology with high-throughput chromosome conformation capture (Hi-C) technique to obtain a chromosome-level assembly of V. radiata cultivar 'KUML4'.

Results

The final assembly of the V. radiata genome was 468.08 Mb in size, with a scaffold N50 of 40.75 Mb. This assembly comprised 11 pseudomolecules, covering 96.94% of the estimated genome size. The genome contained 253.85 Mb (54.76%) of repetitive sequences and 27,667 protein-coding genes. Our gene prediction recovered 98.3% of the highly conserved orthologs based on Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. Comparative analyses using sequence data from single-copy orthologous genes indicated that V. radiata diverged from V. mungo approximately 4.17 million years ago. Moreover, gene family analysis revealed that major gene families associated with defense responses were significantly expanded in V. radiata.

Conclusion

Our chromosome-scale genome assembly of V. radiata cultivar KUML4 will provide a valuable genomic resource, supporting genetic improvement and molecular breeding. This data will also be valuable for future comparative genomics studies among legume species.

Chromosome-scale assembly of the African yam bean genome

Genomics-informed breeding of locally adapted, nutritious, albeit underutilised African crops can help mitigate food and nutrition insecurity challenges in Africa, particularly against the backdrop of climate change. However, utilisation of modern genome-assisted crop improvement tools including genomic selection and genome editing for many African indigenous crops is hampered by the scarcity of genomic resources. Here we report on the assembly of the genome of African yam bean (Sphenostylis stenocarpa), a tuberous legume crop that is indigenous to Africa. By combining Nanopore-based assembly with Hi-C scaffolding, we produced a high-quality chromosome-scale assembly with an N50 of 69.5 Mbp. Using transcriptome evidence from Nanopore RNASeq and protein homology evidence from related crops, we predicted and annotated 31,614 putative protein coding genes. We also show how this genome substantially improves anchoring of genetic markers from African yam bean, confirming its significance as a resource for genetic research in African yam bean.

k-mer genome-wide association study for anthracnose and BCMV resistance in a Phaseolus vulgaris Andean Diversity Panel

Access to broad genomic resources and closely linked marker-trait associations for common beans (Phaseolus vulgaris L.) can facilitate development of improved varieties with increased yield, improved market quality traits, and enhanced disease resistance. The emergence of virulent races of anthracnose (caused by Colletotrichum lindemuthianum) and bean common mosaic virus (BCMV) highlight the need for improved methods to identify and incorporate pan-genomic variation in breeding for disease resistance. We sequenced the P. vulgaris Andean Diversity Panel (ADP) and performed a genome-wide association study (GWAS) to identify associations for resistance to BCMV and eight races of anthracnose. Historical single nucleotide polymorphism (SNP)-chip and phenotypic data enabled a three-way comparison between SNP-chip, reference-based whole genome shotgun sequence (WGS)-SNP, and reference-free k-mer (short nucleotide subsequence) GWAS. Across all traits, there was excellent concordance between SNP-chip, WGS-SNP, and k-mer GWAS results-albeit at a much higher marker resolution for the WGS data sets. Significant k-mer haplotype variation revealed selection of the linked I-gene and Co-u traits in North American breeding lines and cultivars. Due to structural variation, only 9.1 to 47.3% of the significantly associated k-mers could be mapped to the reference genome. Thus, to determine the genetic context of cis-associated k-mers, we generated draft whole genome assemblies of four ADP accessions and identified an expanded local repertoire of disease resistance genes associated with resistance to anthracnose and BCMV. With access to variant data in the context of a pan-genome, high resolution mapping of agronomic traits for common bean is now feasible.

A single-nucleotide polymorphism in PvPW1 encoding β-1,3-glucanase 9 is associated with pod width in Phaseolus vulgaris L

Pod width influences pod size, shape, yield, and consumer preference in snap beans (Phaseolus vulgaris L.). In this study, we map PvPW1, a quantitative trait locus associated with pod width in snap beans, through genotyping and phenotyping of recombinant plants. We identify Phvul.006G072800, encoding the β-1,3-glucanase 9 protein, as the causal gene for PvPW1. The PvPW1G3555 allele is found to positively regulate pod width, as revealed by an association analysis between pod width phenotype and the PvPW1G3555C genotype across 17 bi-parental F2 populations. In total, 97.7% of the 133 wide pod accessions carry PvPW1G3555, while 82.1% of the 78 narrow pod accessions carry PvPW1C3555, indicating strong selection pressure on PvPW1 during common bean breeding. Re-sequencing data from 59 common bean cultivars identify an 8-bp deletion in the intron linked to PvPW1C3555, leading to the development of the InDel marker of PvM436. Genotyping 317 common bean accessions with PvM436 demonstrated that accessions with PvM436247 and PvM436227 alleles have wider pods compared to those with PvM436219 allele, establishing PvM436 as a reliable marker for molecular breeding in snap beans. These findings highlight PvPW1 as a critical gene regulating pod width and underscore the utility of PvM436 in marker-assisted selection for snap bean breeding.

The genome of Lupinus mutabilis: Evolution and genetics of an emerging bio-based crop

Lupinus mutabilis is an under-domesticated legume species from the Andean region of South America. It belongs to the New World lupins clade, which groups several lupin species displaying large genetic variation and adaptability to highly different environments. L. mutabilis is attracting interest as a potential multipurpose crop to diversify the European supply of plant proteins, increase agricultural biodiversity, and fulfill bio-based applications. This study reports the first high-quality L. mutabilis genome assembly, which is also the first sequenced assembly of a New World lupin species. Through comparative genomics and phylogenetics, the evolution of L. mutabilis within legumes and lupins is described, highlighting both genomic similarities and patterns specific to L. mutabilis, potentially linked to environmental adaptations. Furthermore, the assembly was used to study the genetics underlying important traits for the establishment of L. mutabilis as a novel crop, including protein and quinolizidine alkaloids contents in seeds, genomic patterns of classic resistance genes, and genomic properties of L. mutabilis mycorrhiza-related genes. These analyses pointed out copy number variation, differential genomic gene contexts, and gene family expansion through tandem duplications as likely important drivers of the genomic diversity observed for these traits between L. mutabilis and other lupins and legumes. Overall, the L. mutabilis genome assembly will be a valuable resource to conduct genetic research and enable genomic-based breeding approaches to turn L. mutabilis into a multipurpose legume crop.

Convergence and divergence of diploid and tetraploid cotton genomes

Polyploidy is an important driving force in speciation and evolution; however, the genomic basis for parallel selection of a particular trait between polyploids and ancestral diploids remains unexplored. Here we construct graph-based pan-genomes for diploid (A2) and allotetraploid (AD1) cotton species, enabled by an assembly of 50 genomes of genetically diverse accessions. We delineate a mosaic genome map of tetraploid cultivars that illustrates genomic contributions from semi-wild forms into modern cultivars. Pan-genome comparisons identify syntenic and hyper-divergent regions of continued variation between diploid and tetraploid cottons, and suggest an ongoing process of sequence evolution potentially linked to the contrasting genome size change in two subgenomes. We highlight 43% of genetic regulatory relationships for gene expression in diploid encompassing sequence divergence after polyploidy, and specifically characterize six underexplored convergent genetic loci contributing to parallel selection of fiber quality. This study offers a framework for pan-genomic dissection of genetic regulatory components underlying parallel selection of desirable traits in organisms.

Comparative Transcriptomic Analysis Reveals Domestication and Improvement Patterns of Broomcorn Millet (Panicum miliaceum L.)

Broomcorn millet (Panicum miliaceum L.) is one of the earliest crops, domesticated nearly 8000 years ago in northern China. It gradually spread across the entire Eurasian continent, as well as to America and Africa, with recent improvement in various reproductive and vegetative traits. To identify the genes that were selected during the domestication and improvement processes, we performed a comparative transcriptome analysis based on wild types, landraces, and improved cultivars of broomcorn millet at both seeding and filling stages. The variations in gene expression patterns between wild types and landraces and between landraces and improved cultivars were further evaluated to explore the molecular mechanisms underlying the domestication and improvement of broomcorn millet. A total of 2155 and 3033 candidate genes involved in domestication and a total of 84 and 180 candidate genes related to improvement were identified at seedling and filling stages of broomcorn millet, respectively. The annotation results suggested that the genes related to metabolites, stress resistance, and plant hormones were widely selected during both domestication and improvement processes, while some genes were exclusively selected in either domestication or improvement stages, with higher selection pressure detected in the domestication process. Furthermore, some domestication- and improvement-related genes involved in stress resistance either lost their functions or reduced their expression levels due to the trade-offs between stress resistance and productivity. This study provided novel genetic materials for further molecular breeding of broomcorn millet varieties with improved agronomic traits.

T and Z, partial seed coat patterning genes in common bean, provide insight into the structure and protein interactions of a plant MBW complex

Flavonoids are secondary metabolites associated with plant seed coat and flower color. These compounds provide health benefits to humans as anti-inflammatory and antioxidant compounds. The expression of the late biosynthetic genes in the flavonoid pathway is controlled by a ternary MBW protein complex consisting of interfacing MYB, beta-helix-loop-helix (bHLH), and WD40 Repeat (WDR) proteins. P, the master regulator gene of the flavonoid expression in common bean (Phaseolus vulgaris L.), was recently determined to encode a bHLH protein. The T and Z genes control the distribution of color in bean seeds and flowers and have historically been considered regulators of the flavonoid gene expression. T and Z candidates were identified using reverse genetics based on genetic mapping, phylogenetic analysis, and mutant analysis. Domain and AlphaFold2 structure analyses determined that T encodes a seven-bladed β-propeller WDR protein, while Z encodes a R2R3 MYB protein. Deletions and SNPs in T and Z mutants, respectively, altered the 3D structure of these proteins. Modeling of the Z MYB/P bHLH/T WDR MBW complex identified interfacing sequence domains and motifs in all three genes that are conserved in dicots. One Z MYB motif is a possible beta-molecular recognition feature (β-MoRF) that only appears in a structured state when Z MYB is modeled as a component of a MBW complex. Complexes containing mutant T and Z proteins changed the interaction of members of the complex in ways that would alter their role in regulating the expression of genes in the flavonoid pathway.

Loss of Lateral suppressor gene is associated with evolution of root nodule symbiosis in Leguminosae

BACKGROUND

Root nodule symbiosis (RNS) is a fascinating evolutionary event. Given that limited genes conferring the evolution of RNS in Leguminosae have been functionally validated, the genetic basis of the evolution of RNS remains largely unknown. Identifying the genes involved in the evolution of RNS will help to reveal the mystery.

RESULTS

Here, we investigate the gene loss event during the evolution of RNS in Leguminosae through phylogenomic and synteny analyses in 48 species including 16 Leguminosae species. We reveal that loss of the Lateral suppressor gene, a member of the GRAS-domain protein family, is associated with the evolution of RNS in Leguminosae. Ectopic expression of the Lateral suppressor (Ls) gene from tomato and its homolog MONOCULM 1 (MOC1) and Os7 from rice in soybean and Medicago truncatula result in almost completely lost nodulation capability. Further investigation shows that Lateral suppressor protein, Ls, MOC1, and Os7 might function through an interaction with NODULATION SIGNALING PATHWAY 2 (NSP2) and CYCLOPS to repress the transcription of NODULE INCEPTION (NIN) to inhibit the nodulation in Leguminosae. Additionally, we find that the cathepsin H (CTSH), a conserved protein, could interact with Lateral suppressor protein, Ls, MOC1, and Os7 and affect the nodulation.

CONCLUSIONS

This study sheds light on uncovering the genetic basis of the evolution of RNS in Leguminosae and suggests that gene loss plays an essential role.

Genome-Wide Association for Morphological and Agronomic Traits in Phaseolus vulgaris L. Accessions

Exploring genetic resources through genomic analyses has emerged as a powerful strategy to develop common bean (Phaseolus vulgaris L.) cultivars that are both productive and well-adapted to various environments. This study aimed to identify genomic regions linked to morpho-agronomic traits in Mesoamerican and Andean common bean accessions and to elucidate the proteins potentially involved in these traits. We evaluated 109 common bean accessions over three agricultural years, focusing on traits including the grain yield (YDSD), 100-seed weight (SW), number of seeds per pod (SDPD), number of pods per plant (PDPL), first pod insertion height (FPIH), plant height (PLHT), days to flowering (DF), and days to maturity (DPM). Using multilocus methods such as mrMLM, FASTmrMLM, FASTmrEMMA, ISIS EM-BLASSO, and pLARmEB, we identified 36 significant SNPs across all chromosomes (Pv01 to Pv11). Validating these SNPs and candidate genes in segregating populations is crucial for developing more productive common bean cultivars through marker-assisted selection.

Cell size differences affect photosynthetic capacity in a Mesoamerican and an Andean genotype of Phaseolus vulgaris L

The efficiency of CO2 flux in the leaf is hindered by a several structural and biochemical barriers which affect the overall net photosynthesis. However, the dearth of information about the genetic control of these features is limiting our ability for genetic manipulation. We performed a comparative analysis between three-week-old plants of a Mesoamerican and an Andean cultivar of Phaseolus vulgaris at variable light and CO2 levels. The Mesoamerican bean had higher photosynthetic rate, maximum rate of rubisco carboxylase activity and maximum rate of photosynthetic electron transport at light saturation conditions than its Andean counterpart. Leaf anatomy comparison between genotypes showed that the Mesoamerican bean had smaller cell sizes than the Andean bean. Smaller epidermal cells in the Mesoamerican bean resulted in higher stomata density and consequently higher stomatal conductance for water vapor and CO2 than in the Andean bean. Likewise, smaller palisade and spongy mesophyll cells in the Mesoamerican than in the Andean bean increased the cell surface area per unit of volume and consequently increased mesophyll conductance. Finally, smaller cells in the Mesoamerican also increased chlorophyll and protein content per unit of leaf area. In summary, we show that different cell sizes controls the overall net photosynthesis and could be used as a target for genetic manipulation to improve photosynthesis.

GWAS of resistance to three bacterial diseases in the Andean common bean diversity panel

Bacterial brown spot (BBS) caused by Pseudomonas syringae pv. syringae (Pss), common bacterial blight (CBB) caused by Xanthomonas axonopodis pv. phaseoli (Xap) and Xanthomonas fuscans subsp. fuscans (Xff), and halo bacterial blight (HBB), caused by Pseudomonas syringae pv. phaseolicola (Psph), are major bacterial diseases that severely affect common bean yields and global food security. Andean-origin dry beans, representing large-seeded market classes, are particularly susceptible. Using 140,325 SNPs, a multi-locus GWAS was conducted on subsets of the Andean diversity panel (ADP) phenotyped for BBS in South Africa, CBB in Puerto Rico, South Africa, and Zambia, and HBB in South Africa, through natural infection, artificial inoculation, or both. Twenty-four QTL associated with resistance were identified: nine for BBS, eight for CBB, and seven for HBB. Four QTL intervals on Pv01, Pv03, Pv05, and Pv08 overlapped with BBS and HBB resistance. A genomic interval on Pv01, near the fin gene, which determines growth habit, was linked to resistance to all three pathogens. Different QTLs were detected for BBS and CBB resistance when phenotyped under natural infection versus artificial inoculation. These results underscore the importance of combining phenotyping methods in multi-GWAS to capture the full genetic spectrum. Previously recognized CBB resistance QTL SAP6 and SU91 and HBB resistance QTL HB4.2, and HB5.1, were observed. Other common (MAF >0.25) and rare (MAF <0.05) resistance QTL were also detected. Overall, these findings enhance the understanding and utilization of bacterial resistance present in ADP for the development of common beans with improved resistance.

Seed color patterns in domesticated common bean are regulated by MYB-bHLH-WD40 transcription factors and temperature

Seed colors and color patterns are critical for the survival of wild plants and the consumer appeal of crops. In common bean, a major global staple, these patterns are also essential in determining market classes, yet the genetic and environmental control of many pigmentation patterns remains unresolved. In this study, we genetically mapped variation for several important seed pattern loci, including T, Bip, phbw, and Z, which co-segregated with candidate genes PvTTG1, PvMYC1, PvTT8, and PvTT2, respectively. Proteins encoded by these genes are predicted to work together in MYB-bHLH-WD40 (MBW) complexes, propagating flavonoid biosynthesis across the seed coat as observed in Arabidopsis. Whole-genome sequencing of 37 accessions identified mutations, including seven unique parallel mutations in T (PvTTG1) and non-synonymous SNPs in highly conserved residues in bipana (PvMYC1) and z (PvTT2). A 612 bp intron deletion in phbw (PvTT8) eliminated motifs conserved since the Papilionoideae origin and corresponded to a 20-fold reduction in transcript abundance. In multi-location field trials of seven varieties with partial seed coat pigmentation patterning, the pigmented seed coat area correlated positively with ambient temperature, with up to 11-fold increases in the pigmented area from the coolest to the warmest environments. In controlled growth chamber conditions, an increase of 4°C was sufficient to cause pigmentation on an average additional 21% of the seed coat area. Our results shed light on key steps of flavonoid biosynthesis in common bean. They will inform breeding efforts for seed coat color/patterning to improve consumer appeal in this nutritious staple crop.

Antioxidant and Anti-Inflammatory Capacities of Three Dry Bean Varieties after Cooking and In Vitro Gastrointestinal Digestion

The present study delved into the chemical composition, antioxidant, and anti-inflammatory properties of three dry edible beans: Black (BL), Great Northern (GN), and Pinto (PN). The beans were soaked, cooked, and subjected to in vitro gastrointestinal (GI) digestion. BL bean exhibited significantly higher gastric (42%) and intestinal (8%) digestion rates. Comparative assessment of soluble GI-digested fractions (<3 kDa) revealed that the GN bean exhibited the highest abundance of dipeptides (P < 0.05). The BL bean fraction displayed a 4-fold increase in tripeptides (P < 0.05). Both BL and PN bean fractions are high in essential free amino acids, flavonols, and derivatives of hydroxybenzoic acid when compared to the GN bean. All the beans exhibited the ability to mitigate TNF-α-induced pro-inflammatory signaling; however, the BL bean fraction was the most effective at lowering AAPH-induced oxidative stress in HT-29 cells, followed by the GN bean (P < 0.05). In contrast, a low antioxidant effect was observed with PN beans.

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

QTL mapping for pod quality and yield traits in snap bean (Phaseolus vulgaris L.)

Pod quality and yield traits in snap bean (Phaseolus vulgaris L.) influence consumer preferences, crop adoption by farmers, and the ability of the product to be commercially competitive locally and globally. The objective of the study was to identify the quantitative trait loci (QTL) for pod quality and yield traits in a snap × dry bean recombinant inbred line (RIL) population. A total of 184 F6 RILs derived from a cross between Vanilla (snap bean) and MCM5001 (dry bean) were grown in three field sites in Kenya and one greenhouse environment in Davis, CA, USA. They were genotyped at 5,951 single nucleotide polymorphisms (SNPs), and composite interval mapping was conducted to identify QTL for 16 pod quality and yield traits, including pod wall fiber, pod string, pod size, and harvest metrics. A combined total of 44 QTL were identified in field and greenhouse trials. The QTL for pod quality were identified on chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, and Pv07, and for pod yield were identified on Pv08. Co-localization of QTL was observed for pod quality and yield traits. Some identified QTL overlapped with previously mapped QTL for pod quality and yield traits, with several others identified as novel. The identified QTL can be used in future marker-assisted selection in snap bean.

Adaptive gene loss in the common bean pan-genome during range expansion and domestication

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.

De novo genome assembly of white clover (Trifolium repens L.) reveals the role of copy number variation in rapid environmental adaptation

BACKGROUND

White clover (Trifolium repens) is a globally important perennial forage legume. This species also serves as an eco-evolutionary model system for studying within-species chemical defense variation; it features a well-studied polymorphism for cyanogenesis (HCN release following tissue damage), with higher frequencies of cyanogenic plants favored in warmer locations worldwide. Using a newly generated haplotype-resolved genome and two other long-read assemblies, we tested the hypothesis that copy number variants (CNVs) at cyanogenesis genes play a role in the ability of white clover to rapidly adapt to local environments. We also examined questions on subgenome evolution in this recently evolved allotetraploid species and on chromosomal rearrangements in the broader IRLC legume clade.

RESULTS

Integration of PacBio HiFi, Omni-C, Illumina, and linkage map data yielded a completely de novo genome assembly for white clover (created without a priori sequence assignment to subgenomes). We find that white clover has undergone extensive transposon diversification since its origin but otherwise shows highly conserved genome organization and composition with its diploid progenitors. Unlike some other clover species, its chromosomal structure is conserved with other IRLC legumes. We further find extensive evidence of CNVs at the major cyanogenesis loci; these contribute to quantitative variation in the cyanogenic phenotype and to local adaptation across wild North American populations.

CONCLUSIONS

This work provides a case study documenting the role of CNVs in local adaptation in a plant species, and it highlights the value of pan-genome data for identifying contributions of structural variants to adaptation in nature.

Evolution of Einkorn wheat centromeres is driven by the mutualistic interplay of two LTR retrotransposons

BACKGROUND

Centromere function is highly conserved across eukaryotes, but the underlying centromeric DNA sequences vary dramatically between species. Centromeres often contain a high proportion of repetitive DNA, such as tandem repeats and/or transposable elements (TEs). Einkorn wheat centromeres lack tandem repeat arrays and are instead composed mostly of the two long terminal repeat (LTR) retrotransposon families RLG_Cereba and RLG_Quinta which specifically insert in centromeres. However, it is poorly understood how these two TE families relate to each other and if and how they contribute to centromere function and evolution.

RESULTS

Based on conservation of diagnostic motifs (LTRs, integrase and primer binding site and polypurine-tract), we propose that RLG_Cereba and RLG_Quinta are a pair of autonomous and non-autonomous partners, in which the autonomous RLG_Cereba contributes all the proteins required for transposition, while the non-autonomous RLG_Quinta contributes GAG protein. Phylogenetic analysis of predicted GAG proteins showed that the RLG_Cereba lineage was present for at least 100 million years in monocotyledon plants. In contrast, RLG_Quinta evolved from RLG_Cereba between 28 and 35 million years ago in the common ancestor of oat and wheat. Interestingly, the integrase of RLG_Cereba is fused to a so-called CR-domain, which is hypothesized to guide the integrase to the functional centromere. Indeed, ChIP-seq data and TE population analysis show only the youngest subfamilies of RLG_Cereba and RLG_Quinta are found in the active centromeres. Importantly, the LTRs of RLG_Quinta and RLG_Cereba are strongly associated with the presence of the centromere-specific CENH3 histone variant. We hypothesize that the LTRs of RLG_Cereba and RLG_Quinta contribute to wheat centromere integrity by phasing and/or placing CENH3 nucleosomes, thus favoring their persistence in the competitive centromere-niche.

CONCLUSION

Our data show that RLG_Cereba cross-mobilizes the non-autonomous RLG_Quinta retrotransposons. New copies of both families are specifically integrated into functional centromeres presumably through direct binding of the integrase CR domain to CENH3 histone variants. The LTRs of newly inserted RLG_Cereba and RLG_Quinta elements, in turn, recruit and/or phase new CENH3 deposition. This mutualistic interplay between the two TE families and the plant host dynamically maintains wheat centromeres.

Chromosome-level reference genome and resequencing of 322 accessions reveal evolution, genomic imprint and key agronomic traits in adzuki bean

Adzuki bean (Vigna angularis) is an important legume crop cultivated in over 30 countries worldwide. We developed a high-quality chromosome-level reference genome of adzuki bean cultivar Jingnong6 by combining PacBio Sequel long-read sequencing with short-read and Hi-C technologies. The assembled genome covers 97.8% of the adzuki bean genome with a contig N50 of approximately 16 Mb and a total of 32 738 protein-coding genes. We also generated a comprehensive genome variation map of adzuki bean by whole-genome resequencing (WGRS) of 322 diverse adzuki beans accessions including both wild and cultivated. Furthermore, we have conducted comparative genomics and a genome-wide association study (GWAS) on key agricultural traits to investigate the evolution and domestication. GWAS identified several candidate genes, including VaCycA3;1, VaHB15, VaANR1 and VaBm, that exhibited significant associations with domestication traits. Furthermore, we conducted functional analyses on the roles of VaANR1 and VaBm in regulating seed coat colour. We provided evidence for the highest genetic diversity of wild adzuki (Vigna angularis var. nipponensis) in China with the presence of the most original wild adzuki bean, and the occurrence of domestication process facilitating transition from wild to cultigen. The present study elucidates the genetic basis of adzuki bean domestication traits and provides crucial genomic resources to support future breeding efforts in adzuki bean.