Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

Fine mapping of a new common bean anthracnose resistance gene (Co-18) to the proximal end of Pv10 in Indian landrace KRC-5

KEY MESSAGE

Mapping and fine mapping of bean anthracnose resistance genes is a continuous process. We report fine mapping of anthracnose resistance gene Co-18 which is the first anthracnose gene mapped to Pv10. The discovery of resistance gene is a major gain in the bean anthracnose pathosystem research. Among the Indian common bean landraces, KRC-5 exhibit high levels of resistance to the bean anthracnose pathogen Colletotrichum lindemuthianum. To precisely map the anthracnose resistance gene, we used a Recombinant Inbred Line (F2:9 RIL) population (KRC-5 × Jawala). The inheritance test revealed that KRC-5 carries a dominant resistance gene temporarily designated as Co-18. We discovered two RAPD markers linked to Co-18 among 287 RAPD markers. These RAPD markers were eventually developed into SCARs (Sc-OPR15 and Sc-OPF6) and flank Co-18 on chromosome Pv10 at a distance of 5.3 and 4.2 cM, respectively. At 4.0-4.1 Mb on Pv10, we detected a SNP (single-nucleotide polymorphism) signal. We synthesized 58 SSRs and 83 InDels from a pool of 135 SSRs and 1134 InDels, respectively. Five SSRs, four InDels, and two SCARs were used to generate the high-density linkage map, which led to the identification of two SSRs (SSR24 and SSR36) that are tightly linked to Co-18. These two SSRs flank the Co-18 to 178 kb genomic region with 13 candidate genes including five NLR (nucleotide-binding and leucine-rich repeat) genes. The closely linked markers SSR24 and SSR36 will be used in cloning and pyramiding of the Co-18 gene with other R genes to develop durable resistant bean varieties.

Functional, transcriptomic, and lipidomic studies of the choC gene encoding a phospholipid methyltransferase in Aspergillus fumigatus

IMPORTANCE

This study explored the phospholipid metabolic pathway in A. fumigatus and its relationship with fungal growth, metabolism, and pathogenicity. ChoC, based on its critical roles in many aspects of the fungus and relatively conserved characteristics in filamentous fungi with low similarity with mammalian ones, can be a novel target of new antifungal drugs.

Phaseolus vulgaris mutants reveal variation in the nuclear genome

Phaseolus vulgaris L. (common bean) is an essential source of proteins in the human diet worldwide. Bean breeding programs to increase genetic diversity based on induced mutagenesis have a long tradition in Bulgaria. Common bean varieties with high productivity, wide environmental adaptability, good nutritional properties, and improved disease resistance have been successfully developed. In this study, we aimed to investigate selected nuclear genome features, such as the genome size, the number and chromosomal distribution of 5S and 35S rDNA loci by using the fluorescence in situ hybridization (FISH), as well as the level of DNA damage in some local Bulgarian accessions and mutants of P. vulgaris. Flow cytometry analyses revealed no significant differences in genome size between analyzed lines except for one of the analyzed mutants, M19. The value of genome size 2C DNA is about 1.37 pg2C -1 for all lines, whereas it is 1.42 pg2C-1 for M19. The chromosome number remains the same (2n=22) for all analyzed lines. Results of FISH analyses showed that the number of 5S rDNA was stable among accessions and mutant lines (four loci), while the number of 35S rDNA loci was shown as highly polymorphic, varying between ten and sixteen, and displaying differences in the size and location of 35S rDNA loci between analyzed genotypes. The cell cycle profile was different for the analyzed genotypes. The results revealed that wide variation in genome organization and size as well as DNA damage characterizes the analyzed genetic resources of the common bean.

Genome-wide analysis of the SMXL gene family in common bean and identification of karrikin-responsive PvSMXL2 as a negative regulator of PEG-induced drought stress

Common bean (Phaseolus vulgaris L.) is a major legume crop worldwide, but its growth and development frequently face challenges due to abiotic stresses, particularly drought. Proper supplement of copper could mitigate the adverse effects of drought, but excessive accumulation of this metal in plants can be harmful. The suppressor of MAX2 1-like (SMXL) gene family, which plays important roles in various plant processes, including stress responses, remains poorly understood in common bean. In this study, we identified nine orthologues of SMXL genes in common bean, which are located on six chromosomes and classified into four subgroups. Basic molecular properties, including theoretical isoelectric point (PI), molecular weight (MW), grand average of hydropathicity (GVIO), gene structure, and conserved motifs were characterized, and numerous cis-elements in promoters were predicted. The expression patterns of PvSMXL genes were found to be distinct under 10% polyethylene glycol (PEG)-induced drought stress and 200 μM Cu treatments. Most PvSMXLs showed reduced expression in response to Cu treatment, whereas nearly half PvSMXLs exhibited inducible expression under drought stress. PvSMXL2, which exhibited a rapid response to karrikin 1 (KAR1), an active form of the plant growth regulators newly found in the smoke of burning plant material, was down-regulated by both PEG-induced drought and Cu stresses. Transient silencing of PvSMXL2 resulted in enhanced drought stress tolerance without conferring Cu tolerance. These findings provide valuable insights into the functions of SMXL genes in common bean under abiotic stress conditions.

Fine-mapping and evolutionary history of R-BPMV, a dominant resistance gene to Bean pod mottle virus in Phaseolus vulgaris L

KEY MESSAGE

R-BPMV is located within a recently expanded TNL cluster in the Phaseolus genus with suppressed recombination and known for resistance to multiple pathogens including potyviruses controlled by the I gene. Bean pod mottle virus (BPMV) is a comovirus that infects common bean and legumes in general. BPMV is distributed throughout the world and is a major threat on soybean, a closely related species of common bean. In common bean, BAT93 was reported to carry the R-BPMV resistance gene conferring resistance to BPMV and linked with the I resistance gene. To fine map R-BPMV, 182 recombinant inbred lines (RILs) derived from the cross BAT93 × JaloEEP558 were genotyped with polymerase chain reaction (PCR)-based markers developed using genome assemblies from G19833 and BAT93, as well as BAT93 BAC clone sequences. Analysis of RILs carrying key recombination events positioned R-BPMV to a target region containing at least 16 TIR-NB-LRR (TNL) sequences in BAT93. Because the I cluster presents a suppression of recombination and a large number of repeated sequences, none of the 16 TNLs could be excluded as R-BPMV candidate gene. The evolutionary history of the TNLs for the I cluster were reconstructed using microsynteny and phylogenetic analyses within the legume family. A single I TNL was present in Medicago truncatula and lost in soybean, mirroring the absence of complete BPMV resistance in soybean. Amplification of TNLs in the I cluster predates the divergence of the Phaseolus species, in agreement with the emergence of R-BPMV before the separation of the common bean wild centers of diversity. This analysis provides PCR-based markers useful in marker-assisted selection (MAS) and laid the foundation for cloning of R-BPMV resistance gene in order to transfer the resistance into soybean.

Comparative analysis of repetitive DNA in dysploid and non-dysploid Phaseolus beans

Structural karyotype changes result from ectopic recombination events frequently associated with repetitive DNA. Although most Phaseolus species present relatively stable karyotypes with 2n = 22 chromosomes, the karyotypes of species of the Leptostachyus group show high rates of structural rearrangements, including a nested chromosome fusion that led to the dysploid chromosome number of the group (2n = 20). We examined the roles of repetitive landscapes in the rearrangements of species of the Leptostachyus group using genome-skimming data to characterize the repeatome in a range of Phaseolus species and compared them to species of that group (P. leptostachyus and P. macvaughii). LTR retrotransposons, especially the Ty3/gypsy lineage Chromovirus, were the most abundant elements in the genomes. Differences in the abundance of Tekay, Retand, and SIRE elements between P. macvaughii and P. leptostachyus were reflected in their total amounts of Ty3/gypsy and Ty1/copia. The satellite DNA fraction was the most divergent among the species, varying both in abundance and distribution, even between P. leptostachyus and P. macvaughii. The rapid turnover of repeats in the Leptostachyus group may be associated with the several rearrangements observed.

High rates of structural rearrangements have shaped the chromosome evolution in dysploid Phaseolus beans

KEY MESSAGE

Karyotypes evolve through numerical and structural chromosome rearrangements. We show that Phaseolus leptostachyus, a wild bean, underwent a rapid genome reshuffling associated with the reduction from 11 to 10 chromosome pairs, but without whole genome duplication, the highest chromosome evolution rate known for plants. Plant karyotypes evolve through structural rearrangements often associated with polyploidy or dysploidy. The genus Phaseolus comprises ~ 90 species, five of them domesticated due to their nutritional relevance. Most of the species have 2n = 22 karyotypes and are highly syntenic, except for three dysploid karyotypes of species from the Leptostachyus group (2n = 20) that have accumulated several rearrangements. Here, we investigated the degrees of structural rearrangements among Leptostachyus and other Phaseolus groups by estimating their chromosomal evolution rates (CER). For this, we combined our oligo-FISH barcode system for beans and chromosome-specific painting probes for chromosomes 2 and 3, with rDNA and a centromeric probe to establish chromosome orthologies and identify structural rearrangements across nine Phaseolus species. We also integrated the detected rearrangements with a phylogenomic approach to estimate the CERs for each Phaseolus lineage. Our data allowed us to identify translocations, inversions, duplications and deletions, mostly in species belonging to the Leptostachyus group. Phaseolus leptostachyus showed the highest CER (12.31 rearrangements/My), a tenfold increase in contrast to the 2n = 22 species analysed. This is the highest rate known yet for plants, making it a model species for investigating the mechanisms behind rapid genome reshuffling in early species diversification.

Whole-genome resequencing of common bean elite breeding lines

The expansion of bean genome technologies has prompted new perspectives on generating resources and knowledge essential to research and implementing biotechnological tools for the practical operations of plant breeding programs. This study aimed to resequence the entire genome (whole genome sequencing-WGS) of 40 bean genotypes selected based on their significance in breeding programs worldwide, with the objective of generating an extensive database for the identification of single nucleotide polymorphisms (SNPs). Over 6 million SNPs were identified, distributed across the 11 bean chromosomes. After quality variant filtering, 420,509 high-quality SNPs were established, with an average of 38,228 SNPs per chromosome. These variants were categorized based on their predicted effects, revealing that the majority exerted a modifier impact on non-coding genome regions (94.68%). Notably, a significant proportion of SNPs occurred in intergenic regions (62.89%) and at least one SNP was identified in 58.63% of the genes annotated in the bean genome. Of particular interest, 7841 SNPs were identified in 85% of the putative plant disease defense-related genes, presenting a valuable resource for crop breeding efforts. These findings provide a foundation for the development of innovative and broadly applicable technologies for the routine selection of superior genotypes in global bean improvement and germplasm characterization programs.

Impact of drought and salt stress on galactinol and raffinose family oligosaccharides in common bean (Phaseolus vulgaris)

Due to climate change, farmers will face more extreme weather conditions and hence will need crops that are better adapted to these challenges. The raffinose family oligosaccharides (RFOs) could play a role in the tolerance of crops towards abiotic stress. To investigate this, we determined for the first time the importance of galactinol and RFOs in the roots and leaves of common bean under drought and salt stress conditions. Initially, the physiological characteristics of common bean under agronomically relevant abiotic stress conditions were investigated by measuring the growth rate, transpiration rate, chlorophyll concentration and membrane stability, allowing to establish relevant sampling points. Subsequently, the differential gene expression profiles of the galactinol and RFO biosynthetic genes and the amount of galactinol and RFO molecules were measured in the primary leaves and roots of Phaseolus vulgaris cv. CIAP7247F at these sampling points, using RT-qPCR and HPAEC-PAD, respectively. Under drought stress, the genes galactinol synthase 1, galactinol synthase 3 and stachyose synthase were significantly upregulated in the leaves and had a high transcript level in comparison with the other galactinol and RFO biosynthetic genes. This was in accordance with the significantly higher amount of galactinol and raffinose detected in the leaves. Under salt stress, raffinose was also present in a significantly higher quantity in the leaves. In the roots, transcript levels of the RFO biosynthetic genes were generally low and no galactinol, raffinose or stachyose could be detected. These results suggest that in the leaves, both galactinol and raffinose could play a role in the protection of common bean against abiotic stresses. Especially, the isoform galactinol synthase 3 could have a specific role during drought stress and forms an interesting candidate to improve the abiotic stress resistance of common bean or other plant species.

A comparative analysis of stably expressed genes across diverse angiosperms exposes flexibility in underlying promoter architecture

Promoters regulate both the amplitude and pattern of gene expression-key factors needed for optimization of many synthetic biology applications. Previous work in Arabidopsis found that promoters that contain a TATA-box element tend to be expressed only under specific conditions or in particular tissues, while promoters which lack any known promoter elements, thus designated as Coreless, tend to be expressed more ubiquitously. To test whether this trend represents a conserved promoter design rule, we identified stably expressed genes across multiple angiosperm species using publicly available RNA-seq data. Comparisons between core promoter architectures and gene expression stability revealed differences in core promoter usage in monocots and eudicots. Furthermore, when tracing the evolution of a given promoter across species, we found that core promoter type was not a strong predictor of expression stability. Our analysis suggests that core promoter types are correlative rather than causative in promoter expression patterns and highlights the challenges in finding or building constitutive promoters that will work across diverse plant species.

Identification and Characterization of Common Bean (Phaseolus vulgaris) Non-Nodulating Mutants Altered in Rhizobial Infection

The symbiotic N₂-fixation process in the legume-rhizobia interaction is relevant for sustainable agriculture. The characterization of symbiotic mutants, mainly in model legumes, has been instrumental for the discovery of symbiotic genes, but similar studies in crop legumes are scant. To isolate and characterize common bean (Phaseolus vulgaris) symbiotic mutants, an ethyl methanesulphonate-induced mutant population from the BAT 93 genotype was analyzed. Our initial screening of Rhizobium etli CE3-inoculated mutant plants revealed different alterations in nodulation. We proceeded with the characterization of three non-nodulating (nnod), apparently monogenic/recessive mutants: nnod(1895), nnod(2353) and nnod(2114). Their reduced growth in a symbiotic condition was restored when the nitrate was added. A similar nnod phenotype was observed upon inoculation with other efficient rhizobia species. A microscopic analysis revealed a different impairment for each mutant in an early symbiotic step. nnod(1895) formed decreased root hair curling but had increased non-effective root hair deformation and no rhizobia infection. nnod(2353) produced normal root hair curling and rhizobia entrapment to form infection chambers, but the development of the latter was blocked. nnod(2114) formed infection threads that did not elongate and thus did not reach the root cortex level; it occasionally formed non-infected pseudo-nodules. The current research is aimed at mapping the responsible mutated gene for a better understanding of SNF in this critical food crop.

Mining legume germplasm for genetic gains: An Indian perspective

Legumes play a significant role in food and nutritional security and contribute to environmental sustainability. Although legumes are highly beneficial crops, it has not yet been possible to enhance their yield and production to a satisfactory level. Amid a rising population and low yield levels, per capita average legume consumption in India has fallen by 71% over the last 50 years, and this has led to protein-related malnutrition in a large segment of the Indian population, especially women and children. Several factors have hindered attempts to achieve yield enhancement in grain legumes, including biotic and abiotic pressures, a lack of good ideotypes, less amenability to mechanization, poorer responsiveness to fertilizer input, and a poor genetic base. Therefore, there is a need to mine the approximately 0.4 million ex situ collections of legumes that are being conserved in gene banks globally for identification of ideal donors for various traits. The Indian National Gene Bank conserves over 63,000 accessions of legumes belonging to 61 species. Recent initiatives have been undertaken in consortia mode with the aim of unlocking the genetic potential of ex situ collections and conducting large-scale germplasm characterization and evaluation analyses. We assume that large-scale phenotyping integrated with omics-based science will aid the identification of target traits and their use to enhance genetic gains. Additionally, in cases where the genetic base of major legumes is narrow, wild relatives have been evaluated, and these are being exploited through pre-breeding. Thus far, >200 accessions of various legumes have been registered as unique donors for various traits of interest.

Comparative Transcriptome Analysis of Tolerant and Sensitive Genotypes of Common Bean (Phaseolus vulgaris L.) in Response to Terminal Drought Stress

We conducted a genome-wide transcriptomic analysis of three drought tolerant and sensitive genotypes of common bean to examine their transcriptional responses to terminal drought stress. We then conducted pairwise comparisons between the root and leaf transcriptomes from the resulting tissue based on combined transcriptomic data from the tolerant and sensitive genotypes. Our transcriptomic data revealed that 491 (6.4%) DEGs (differentially expressed genes) were upregulated in tolerant genotypes, whereas they were downregulated in sensitive genotypes; likewise, 396 (5.1%) DEGs upregulated in sensitive genotypes were downregulated in tolerant genotypes. Several transcription factors, heat shock proteins, and chaperones were identified in the study. Several DEGs in drought DB (data Base) overlapped between genotypes. The GO (gene ontology) terms for biological processes showed upregulation of DEGs in tolerant genotypes for sulfate and drug transmembrane transport when compared to sensitive genotypes. A GO term for cellular components enriched with upregulated DEGs for the apoplast in tolerant genotypes. These results substantiated the temporal pattern of root growth (elongation and initiation of root growth), and ABA-mediated drought response in tolerant genotypes. KEGG (kyoto encyclopedia of genes and genomes) analysis revealed an upregulation of MAPK (mitogen activated protein kinase) signaling pathways and plant hormone signaling pathways in tolerant genotypes. As a result of this study, it will be possible to uncover the molecular mechanisms of drought tolerance in response to terminal drought stress in the field. Further, genome-wide transcriptomic analysis of both tolerant and sensitive genotypes will assist us in identifying potential genes that may contribute to improving drought tolerance in the common bean.

Exchanges of economic plants along the land silk road

BACKGROUD

The greatest contribution of the Silk Road is to communicate among different countries and nationalities, and promote two-way cultural exchanges between the East and the West. We now have clearer understanding about how material civilization and religious culture of Central Asia and West Asia spread eastward along the Land Silk Road. However, there is controversial about how crops migrate along the Land Silk Road.

RESULTS

We summarize archaeology, genetics, and genomics data to explore crop migration patterns. Of the 207 crops that were domesticated along the Land Silk Road, 19 for which genomic evidence was available were selected for discussion.

CONCLUSIONS

There were conflicting lines of evidence for the domestication of Tibetan barley, mustard, lettuce, buckwheat, and chickpea. The main reasons for the conflicting results may include incomplete early knowledge, record differences in different period, sample sizes, and data analysis techniques. There was strong evidence that Tibetan barley, barley, wheat, and jujube were introduced into China before the existence of the Land Silk Road; and mustard, lettuce, buckwheat, chickpea, alfalfa, walnut, cauliflower, grape, spinach, apple, cucumber, mulberry, and pea spread to China via trade and human migration along the Land Silk Road.

Advancing Grain Legumes Domestication and Evolution Studies with Genomics

Grain legumes were domesticated in parallel with cereals in several regions of the world and formed the economic basis of early farming cultures. Since then, legumes have played a vital role in human and animal diets and in fostering agrobiodiversity. Increasing grain legume cultivation will be crucial to safeguard nutritional security and the resilience of agricultural ecosystems across the globe. A better understanding of the molecular underpinnings of domestication and crop evolution of grain legumes may be translated into practical approaches in modern breeding programs to stabilize yield, which is threatened by evolving pathogens and changing climates. During recent decades, domestication research in all crops has greatly benefited from the fast progress in genomic technologies. Yet still, many questions surrounding the domestication and diversification of legumes remain unanswered. In this review, we assess the potential of genomic approaches in grain legume research. We describe the centers of origin and the crucial domestication traits of grain legumes. In addition, we survey the effect of domestication on both above-ground and below-ground traits that have economic importance. Finally, we discuss open questions in grain legume domestication and diversification and outline how to bridge the gap between the preservation of historic crop diversity and their utilization in modern plant breeding.

Escalate protein plates from legumes for sustainable human nutrition

Protein is one of the most important, foremost, and versatile nutrients in food. The quantity and quality of protein are determinants of its nutritional values. Therefore, adequate consumption of high-quality protein is essential for optimal growth, development, and health of humans. Based on short-term nitrogen balance studies, the Recommended Dietary Allowance of protein for the healthy adult with minimal physical activity is 0.8 g protein/kg body weight (BW) per day. Proteins are present in good quantities in not only animals but also in plants, especially in legumes. With the growing demand for protein, interest in plant proteins is also rising due to their comparative low cost as well as the increase in consumers' demand originating from health and environmental concerns. Legumes are nutrient-dense foods, comprising components identified as "antinutritional factors" that can reduce the bioavailability of macro and micronutrients. Other than nutritive value, the physiochemical and behavioral properties of proteins during processing plays a significant role in determining the end quality of food. The term "complete protein" refers to when all nine essential amino acids are present in the correct proportion in our bodies. To have a balanced diet, the right percentage of protein is required for our body. The consumption of these high protein-containing foods will lead to protein sustainability and eradicate malnutrition. Here, we shed light on major opportunities to strengthen the contribution of diversity in legume crops products to sustainable diets. This review will boost awareness and knowledge on underutilized proteinous foods into national nutritional security programs.

Characterization of fungal pathogens and germplasm screening for disease resistance in the main production area of the common bean in Argentina

The common bean (Phaseolus vulgaris L.) is the most important grain legume in the human diet, mainly in Africa and Latin America. Argentina is one of the five major producers of the common bean in the world, and the main cultivation areas are concentrated in the northwestern provinces of this country. Crop production of the common bean is often affected by biotic factors like some endemic fungal diseases, which exert a major economic impact on the region. The most important fungal diseases affecting the common bean in Argentina are white mold caused by Sclerotinia sclerotiorum, angular leaf spot caused by Pseudocercospora griseola, web blight and root rot caused by Rhizoctonia solani, which can cause production losses of up to 100% in the region. At the present, the most effective strategy for controlling these diseases is the use of genetic resistance. In this sense, population study and characterization of fungal pathogens are essential for developing cultivars with durable resistance. In this review we report diversity studies carried out on these three fungal pathogens affecting the common bean in northwestern Argentina, analyzing more than 200 isolates by means of molecular, morphological and pathogenic approaches. Also, the screening of physiological resistance in several common bean commercial lines and wild native germplasm is reviewed. This review contributes to the development of sustainable management strategies and cultural practices in bean production aimed to minimize yield losses due to fungal diseases in the common bean.

Genome-wide identification, characterization and expression analysis of glutaredoxin gene family (Grxs) in Phaseolus vulgaris

Glutaredoxins (Grxs) are ubiquitous oxidoreductase proteins implicated in development and abiotic stress response mainly through maintaining redox homoeostasis. Here, we conducted the first systematic analysis of the Grx gene family (PvGrx) in the most popular legume Phaseolus vulgaris (common bean). A total of 50 PvGrx genes were identified, and divided into four classes (CC-type, CGFS-type, CPYC-type and Grl-type) based on the phylogenetic analysis. The different classes have different introns-exons structures and conserved motifs, indicating functional divergence in the PvGrx family. Both tandem and segmental duplications were found to be involved in the expansion of PvGrx family that underwent a purifying selection by excluding the deleterious loss-of-function mutations. Cis-acting regulatory elements and gene ontology analyses predicted their role of distinctive members in abiotic stress response and hormonal signalling. RNA-seq based expression analysis revealed their differential expression pattern during plant development. On the other hand, RT q-PCR analysis revealed that target PvGrx isoforms were associated with nodule organogenesis and symbiosis based on their expression profiles. In addition, a battery of PvGrx candidates were markedly upregulated by different abiotic stressors suggesting their broad spectrum of functions. These findings serve as a reference for functional analysis and genetic improvement in P. vulgaris and related legume species.

Genetic Augmentation of Legume Crops Using Genomic Resources and Genotyping Platforms for Nutritional Food Security

Recent advances in next generation sequencing (NGS) technologies have led the surge of genomic resources for the improvement legume crops. Advances in high throughput genotyping (HTG) and high throughput phenotyping (HTP) enable legume breeders to improve legume crops more precisely and efficiently. Now, the legume breeder can reshuffle the natural gene combinations of their choice to enhance the genetic potential of crops. These genomic resources are efficiently deployed through molecular breeding approaches for genetic augmentation of important legume crops, such as chickpea, cowpea, pigeonpea, groundnut, common bean, lentil, pea, as well as other underutilized legume crops. In the future, advances in NGS, HTG, and HTP technologies will help in the identification and assembly of superior haplotypes to tailor the legume crop varieties through haplotype-based breeding. This review article focuses on the recent development of genomic resource databases and their deployment in legume molecular breeding programmes to secure global food security.

Genome-Wide Association Study and Genomic Prediction for Bacterial Wilt Resistance in Common Bean (Phaseolus vulgaris) Core Collection

Common bean (Phaseolus vulgaris) is one of the major legume crops cultivated worldwide. Bacterial wilt (BW) of common bean (Curtobacterium flaccumfaciens pv. flaccumfaciens), being a seed-borne disease, has been a challenge in common bean producing regions. A genome-wide association study (GWAS) was conducted to identify SNP markers associated with BW resistance in the USDA common bean core collection. A total of 168 accessions were evaluated for resistance against three different isolates of BW. Our study identified a total of 14 single nucleotide polymorphism (SNP) markers associated with the resistance to BW isolates 528, 557, and 597 using mixed linear models (MLMs) in BLINK, FarmCPU, GAPIT, and TASSEL 5. These SNPs were located on chromosomes Phaseolus vulgaris [Pv]02, Pv04, Pv08, and Pv09 for isolate 528; Pv07, Pv10, and Pv11 for isolate 557; and Pv04, Pv08, and Pv10 for isolate 597. The genomic prediction accuracy was assessed by utilizing seven GP models with 1) all the 4,568 SNPs and 2) the 14 SNP markers. The overall prediction accuracy (PA) ranged from 0.30 to 0.56 for resistance against the three BW isolates. A total of 14 candidate genes were discovered for BW resistance located on chromosomes Pv02, Pv04, Pv07, Pv08, and Pv09. This study revealed vital information for developing genetic resistance against the BW pathogen in common bean. Accordingly, the identified SNP markers and candidate genes can be utilized in common bean molecular breeding programs to develop novel resistant cultivars.

The genomic signature of wild‐to‐crop introgression during the domestication of scarlet runner bean ( Phaseolus coccineus L.)

The scarlet runner bean (Phaseolus coccineus) is one of the five domesticated Phaseolus species. It is cultivated in small‐scale agriculture in the highlands of Mesoamerica for its dry seeds and immature pods, and unlike the other domesticated beans, P. coccineus is an open‐pollinated legume. Contrasting with its close relative, the common bean, few studies focusing on its domestication history have been conducted. Demographic bottlenecks associated with domestication might reduce genetic diversity and facilitate the accumulation of deleterious mutations. Conversely, introgression from wild relatives could be a source of variation. Using Genotyping by Sequencing data (79,286 single‐nucleotide variants) from 237 cultivated and wild samples, we evaluated the demographic history of traditional varieties from different regions of Mexico and looked for evidence of introgression between sympatric wild and cultivated populations. Traditional varieties have high levels of diversity, even though there is evidence of a severe initial genetic bottleneck followed by a population expansion. Introgression from wild to domesticated populations was detected, which might contribute to the recovery of the genetic variation. Introgression has occurred at different times: constantly in the center of Mexico; recently in the North West; and anciently in the South. Several factors are acting together to increase and maintain genetic diversity in P. coccineus cultivars, such as demographic expansion and introgression. Wild relatives represent a valuable genetic resource and have played a key role in scarlet runner bean evolution via introgression into traditional varieties.

Genome-Wide Analysis of Soybean Lateral Organ Boundaries Domain Gene Family Reveals the Role in Phytophthora Root and Stem Rot

The plant-specific lateral organ boundaries (LOB) domain (LBD) proteins, a family of transcription factors, play important roles in plant growth and development, as well as in responses to various stresses. However, little is known about the functions of LBD genes in soybean (Glycine max). In this study, we investigated the evolution and classification of the LBD family in soybean by a phylogenetic tree of the LBD gene family from 16 species. Phylogenetic analysis categorized these proteins into two classes (Class I and Class II) with seven subgroups. Moreover, we found that all the 18 LBD ancestors in angiosperm were kept in soybean, common bean genomes, and genome-wide duplication, suggesting the main force for the expansion of LBD from common bean to soybean. Analysis of gene expression profiling data indicated that 16 GmLBD genes were significantly induced at different time points after inoculation of soybean plants (cv. Huachun 6) with Phytophthora sojae (P. sojae). We further assessed the role of four highly upregulated genes, GmLBD9, GmLBD16, GmLBD23, and GmLBD88, in plant defense in soybean hairy roots using the transient overexpression and knockdown assays. The results showed that GmLBD9 and GmLBD23 negatively regulate plant immunity against P. sojae, whereas GmLBD16 and GmLBD88 positively manipulate plant immunity against P. sojae. Collectively, our findings expand our knowledge of the origin and evolution of the GmLBD gene family in soybean and promote the potential application of these genes in soybean genetic improvement.

Identification of conserved and new miRNAs that affect nodulation and strain selectivity in the Phaseolus vulgaris-Rhizobium etli symbiosis through differential analysis of host small RNAs

Phaseolus vulgaris plants from the Mesoamerican centre of genetic diversification establish a preferential and more efficient root nodule symbiosis with sympatric Rhizobium etli strains. This is mediated by changes in host gene expression, which might occur either at the transcriptional or at the post-transcriptional level. However, the implication of small RNA (sRNA)-mediated control of gene expression in strain selectivity has remained elusive. sRNA sequencing was used to identify host microRNAs (miRNAs) differentially regulated in roots at an early stage of the symbiotic interaction, which were further characterized by applying a reverse genetic approach. In silico analysis identified known and new miRNAs that accumulated to a greater extent in the preferential and more efficient interaction. One of them, designated as Pvu-miR5924, participates in the mechanisms that determine the selection of R. etli strains that will colonize the nodules. In addition, the functional analysis of Pvu-miR390b verified that this miRNA is a negative modulator of nodule formation and bacterial infection. This study not only extended the list of miRNAs identified in P. vulgaris but also enabled the identification of miRNAs that play relevant functions in nodule formation, rhizobial infection and the selection of the rhizobial strains that will occupy the nodule.

Genetic Control of Seed Phytate Accumulation and the Development of Low-Phytate Crops: A Review and Perspective

Breeding low phytic acid (lpa) crops is a strategy that has potential to both improve the nutritional quality of food and feed and contribute to the sustainability of agriculture. Here, we review the lipid-independent and -dependent pathways of phytate synthesis and their regulatory mechanisms in plants. We compare the genetic variation of the phytate concentration and distribution in seeds between dicot and monocot species as well as the associated temporal and spatial expression patterns of the genes involved in phytate synthesis and transport. Quantitative trait loci or significant single nucleotide polymorphisms for the seed phytate concentration have been identified in different plant species by linkage and association mapping, and some genes have been cloned from lpa mutants. We summarize the effects of various lpa mutations on important agronomic traits in crop plants and propose SULTR3;3 and SULTR3;4 as optimal target genes for lpa crop breeding.

Nodulation competitiveness and diversification of symbiosis genes in common beans from the American centers of domestication

Phaseolus vulgaris (common bean), having a proposed Mexican origin within the Americas, comprises three centers of diversification: Mesoamerica, the southern Andes, and the Amotape-Huancabamba Depression in Peru-Ecuador. Rhizobium etli is the predominant rhizobium found symbiotically associated with beans in the Americasalthough closely related Rhizobium phylotypes have also been detected. To investigate if symbiosis between bean varieties and rhizobia evolved affinity, firstly nodulation competitiveness was studied after inoculation with a mixture of sympatric and allopatric rhizobial strains isolated from the respective geographical regions. Rhizobia strains harboring nodC types α and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upgamma$$\end{document}γ, which were found predominant in Mexico and Ecuador, were comparable in nodule occupancy at 50% of each in beans from the Mesoamerican and Andean gene pools, but it is one of those two nodC types which clearly predominated in Ecuadorian-Peruvian beans as well as in Andean beans nodC type \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upgamma$$\end{document}γ predominated the sympatric nodC type δ. The results indicated that those beans from Ecuador-Peru and Andean region, respectively exhibited no affinity for nodulation by the sympatric rhizobial lineages that were found to be predominant in bean nodules formed in those respective areas. Unlike the strains isolated from Ecuador, Rhizobium etli isolated from Mexico as well from the southern Andes was highly competitive for nodulation in beans from Ecuador-Peru, and quite similarly competitive in Mesoamerican and Andean beans. Finally, five gene products associated with symbiosis were examined to analyze variations that could be correlated with nodulation competitiveness. A small GTPase RabA2, transcriptional factors NIN and ASTRAY, and nodulation factor receptors NFR1 and NFR5- indicated high conservation but NIN, NFR1 and NFR5 of beans representative of the Ecuador-Peru genetic pool clustered separated from the Mesoamerican and Andean showing diversification and possible different interaction. These results indicated that both host and bacterial genetics are important for mutual affinity, and that symbiosis is another trait of legumes that could be sensitive to evolutionary influences and local adaptation.

Transcriptomic analysis of Mesoamerican and Andean Phaseolus vulgaris accessions revealed mRNAs and lncRNAs associated with strain selectivity during symbiosis

Legume plants establish a nitrogen-fixing symbiosis with soil bacteria known as rhizobia. Compatibility between legumes and rhizobia is determined at species-specific level, but variations in the outcome of the symbiotic process are also influenced by the capacity of the plant to discriminate and select specific strains that are better partners. We compared the transcriptional response of two genetically diverse accessions of Phaseolus vulgaris from Mesoamerica and South Andes to Rhizobium etli strains that exhibit variable degrees of symbiotic affinities. Our results indicate that the plant genotype is the major determinant of the transcriptional reprogramming occurring in roots at early stages of the symbiotic interaction. Differentially expressed genes (DEGs) regulated in the Mesoamerican and the Andean accessions in response to specific strains are different, but they belong to the same functional categories. The common and strain-specific transcriptional responses to rhizobia involve distinct transcription factors  and cis-elements present in the promoters of DEGs in each accession, showing that diversification and domestication of common bean at different geographic regions influenced the evolution of symbiosis differently in each genetic pool. Quantitative PCR analysis validated our transcriptional datasets, which constitute a valuable source of coding and non-coding candidate genes to further unravel the molecular determinants governing the mechanisms by which plants select bacterial strains that produce a better symbiotic outcome.

Phaseolus vulgaris-Colletotrichum lindemuthianum Pathosystem in the Post-Genomic Era: An Update

Phaseolus vulgaris-Colletotrichum lindemuthianum is one among the oldest host and pathogen interface. Researchers have taken painstaking efforts across the world for understanding the dialogue during early and late phases of interaction. Collectively, these efforts resulted in the deluge of information that helped the researchers to underpin the interface. The latest molecular biology techniques furnished novel detection methods for the anthracnose pathogen, refined the understanding of pathogen population dynamics, and provided the insights on co-evolutionary common bean resistance and C. lindemuthianum virulence dynamics. One of the important breakthroughs came when the Phaseolus vulgaris and its corresponding anthracnose pathogen (C. lindemuthianum) genomes were decoded in 2014 and 2017, respectively. Availability of both the genomes yielded a significant genomic information that helped bean communities to fine map the economically important traits and to identify the pathogenicity determinants and effector molecules. The interface is in a continuous development as knowledge of the anthracnose resistance genes, their precise physical locations, and the identification of effector proteins; the fungus arsenals are being routinely updated. Hence, we revisited the interface and tried to provide an overview of host pathogen dialogue in the genomic era. Additionally, we compiled the sporadic information on this pathosystem from India and provided its futuristic road map to shape its research in the world and northern India, the major dry bean area in the country.

Gramene: A Resource for Comparative Analysis of Plants Genomes and Pathways

Gramene is an integrated bioinformatics resource for accessing, visualizing, and comparing plant genomes and biological pathways. Originally targeting grasses, Gramene has grown to host annotations for over 90 plant genomes including agronomically important cereals (e.g., maize, sorghum, wheat, teff), fruits and vegetables (e.g., apple, watermelon, clementine, tomato, cassava), specialty crops (e.g., coffee, olive tree, pistachio, almond), and plants of special or emerging interest (e.g., cotton, tobacco, cannabis, or hemp). For some species, the resource includes multiple varieties of the same species, which has paved the road for the creation of species-specific pan-genome browsers. The resource also features plant research models, including Arabidopsis and C4 warm-season grasses and brassicas, as well as other species that fill phylogenetic gaps for plant evolution studies. Its strength derives from the application of a phylogenetic framework for genome comparison and the use of ontologies to integrate structural and functional annotation data. This chapter outlines system requirements for end-users and database hosting, data types and basic navigation within Gramene, and provides examples of how to (1) explore Gramene's search results, (2) explore gene-centric comparative genomics data visualizations in Gramene, and (3) explore genetic variation associated with a gene locus. This is the first publication describing in detail Gramene's integrated search interface-intended to provide a simplified entry portal for the resource's main data categories (genomic location, phylogeny, gene expression, pathways, and external references) to the most complete and up-to-date set of plant genome and pathway annotations.

Integrated multi-omics analysis provides insights into genome evolution and phosphorus deficiency adaptation in pigeonpea (Cajanus cajan)

Pigeonpea (Cajanus cajan) is an important legume food crop and plays a crucial role in a secure food supply in many developing countries. Several previous studies have suggested that pigeonpea has great potential for phosphorus (P) deficiency tolerance, but little is known about the underlying mechanism. In this study, the physiological and molecular responses of pigeonpea roots to phosphate (Pi) starvation were investigated through integrating phenotypic, genomic, transcriptomic, metabolomic, and lipidomic analyses. The results showed that low-Pi treatment increased total root length, root surface area, and root acid phosphatase activity, and promoted the secretion of organic acids (e.g. citric acids, piscidic acids, and protocatechuic acids) and the degradation of phospholipids and other P-containing metabolites in the roots of pigeonpea. Consistent with the morphological, physiological, and biochemical changes, a large number of genes involved in these Pi-starvation responses were significantly upregulated in Pi-deficient pigeonpea roots. Among these Pi-starvation response genes upregulated by low-Pi treatment, four gene families were expanded through recent tandem duplication in the pigeonpea genome, namely phosphate transporter 1 (PHT1), phosphoethanolamine/phosphocholine phosphatase (PECP), fasciclin-like arabinogalactan protein (FLA), and glutamate decarboxylase (GAD). These gene families may be associated with Pi uptake from the soil, phospholipid recycling, root morphological remodeling, and regulation of organic acid exudation. Taken together, our results suggest that pigeonpea employs complex Pi-starvation responses to strengthen P acquisition and utilization during low-Pi stress. This study provides new insights into the genome evolution and P deficiency adaptation mechanism of pigeonpea.

The molecular events underpinning cultivar differences in melatonin counteracting salt damage in Phaseolus vulgaris

Melatonin (N -acetyl-5-methoxytryptamine) plays important roles in multiple stress responses, especially under salt stress. However, cultivar differences in melatonin mediated salt stress tolerance are unclear. Phaseolus vulgaris L. (common bean) cultivars Jiyin 1 (JY, salt-tolerant) and Xuliyabai (XL, salt-sensitive) were used in this study. Exogenous melatonin significantly improved root growth under salt stress in JY, but had little effect on XL. Physiology analysis showed significant differences in activities of antioxidant enzymes (superoxide, SOD; and catalase, CAT) and malondialdehyde content between JY and XL. Meanwhile, the change of ABA content in JY and XL root was opposite in salt plus melatonin treatment. Comparative root transcriptomes of JY and XL revealed 3505 and 668 differentially expressed genes (DEGs) regulated by salt stress and melatonin. The most enriched melatonin-responsive genes under salt stress are mainly involved in regulation of transcription, oxidation-reduction process, transcription factor activity, oxidoreductase activity. In addition, melatonin induced more obvious changes of DEGs in JY than that in XL under salt condition. Melatonin also significantly induced 41 DEGs only in JY, including signal transduction genes, transcription factors, ubiquitin protein ligases, ion homeostasis and osmotic adjustment genes etc. This study uncovered the molecular mechanism of cultivar difference of melatonin response under salt stress in common bean.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

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

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

Common bacterial blight of bean: a model of seed transmission and pathological convergence

BACKGROUND

Xanthomonas citri pv. fuscans (Xcf) and Xanthomonas phaseoli pv. phaseoli (Xpp) are the causal agents of common bacterial blight of bean (CBB), an important disease worldwide that remains difficult to control. These pathogens belong to distinct species within the Xanthomonas genus and have undergone a dynamic evolutionary history including the horizontal transfer of genes encoding factors probably involved in adaptation to and pathogenicity on common bean. Seed transmission is a key point of the CBB disease cycle, favouring both vertical transmission of the pathogen and worldwide distribution of the disease through global seed trade.

TAXONOMY

Kingdom: Bacteria; phylum: Proteobacteria; class: Gammaproteobacteria; order: Lysobacterales (also known as Xanthomonadales); family: Lysobacteraceae (also known as Xanthomonadaceae); genus: Xanthomonas; species: X. citri pv. fuscans and X. phaseoli pv. phaseoli (Xcf-Xpp).

HOST RANGE

The main host of Xcf-Xpp is the common bean (Phaseolus vulgaris). Lima bean (Phaseolus lunatus) and members of the Vigna genus (Vigna aconitifolia, Vigna angularis, Vigna mungo, Vigna radiata, and Vigna umbellata) are also natural hosts of Xcf-Xpp. Natural occurrence of Xcf-Xpp has been reported for a handful of other legumes such as Calopogonium sp., Pueraria sp., pea (Pisum sativum), Lablab purpureus, Macroptilium lathyroides, and Strophostyles helvola. There are conflicting reports concerning the natural occurrence of CBB agents on tepary bean (Phaseolus acutifolius) and cowpea (Vigna unguiculata subsp. unguiculata).

SYMPTOMS

CBB symptoms occur on all aerial parts of beans, that is, seedlings, leaves, stems, pods, and seeds. Symptoms initially appear as water-soaked spots evolving into necrosis on leaves, pustules on pods, and cankers on twigs. In severe infections, defoliation and wilting may occur.

DISTRIBUTION

CBB is distributed worldwide, meaning that it is frequently encountered in most places where bean is cultivated in the Americas, Asia, Africa, and Oceania, except for arid tropical areas. Xcf-Xpp are regulated nonquarantine pathogens in Europe and are listed in the A2 list by the European and Mediterranean Plant Protection Organization (EPPO).

GENOME

The genome consists of a single circular chromosome plus one to four extrachromosomal plasmids of various sizes, for a total mean size of 5.27 Mb with 64.7% GC content and an average predicted number of 4,181 coding sequences.

DISEASE CONTROL

Management of CBB is based on integrated approaches that comprise measures aimed at avoiding Xcf-Xpp introduction through infected seeds, cultural practices to limit Xcf-Xpp survival between host crops, whenever possible the use of tolerant or resistant bean genotypes, and chemical treatments, mainly restricted to copper compounds. The use of pathogen-free seeds is essential in an effective management strategy and requires appropriate sampling, detection, and identification methods. USEFUL WEBSITES: https://gd.eppo.int/taxon/XANTPH, https://gd.eppo.int/taxon/XANTFF, and http://www.cost.eu/COST_Actions/ca/CA16107.

Chromosome-length genome assemblies of six legume species provide insights into genome organization, evolution, and agronomic traits for crop improvement

INTRODUCTION

Legume crops are an important source of protein and oil for human health and in fixing atmospheric N₂ for soil enrichment. With an objective to accelerate much-needed genetic analyses and breeding applications, draft genome assemblies were generated in several legume crops; many of them are not high quality because they are mainly based on short reads. However, the superior quality of genome assembly is crucial for a detailed understanding of genomic architecture, genome evolution, and crop improvement.

OBJECTIVES

Present study was undertaken with an objective of developing improved chromosome-length genome assemblies in six different legumes followed by their systematic investigation to unravel different aspects of genome organization and legume evolution.

METHODS

We employed in situ Hi-C data to improve the existing draft genomes and performed different evolutionary and comparative analyses using improved genome assemblies.

RESULTS

We have developed chromosome-length genome assemblies in chickpea, pigeonpea, soybean, subterranean clover, and two wild progenitor species of cultivated groundnut (A. duranensis and A. ipaensis). A comprehensive comparative analysis of these genome assemblies offered improved insights into various evolutionary events that shaped the present-day legume species. We highlighted the expansion of gene families contributing to unique traits such as nodulation in legumes, gravitropism in groundnut, and oil biosynthesis in oilseed legume crops such as groundnut and soybean. As examples, we have demonstrated the utility of improved genome assemblies for enhancing the resolution of "QTL-hotspot" identification for drought tolerance in chickpea and marker-trait associations for agronomic traits in pigeonpea through genome-wide association study. Genomic resources developed in this study are publicly available through an online repository, 'Legumepedia'.

CONCLUSION

This study reports chromosome-length genome assemblies of six legume species and demonstrates the utility of these assemblies in crop improvement. The genomic resources developed here will have significant role in accelerating genetic improvement applications of legume crops.

Quinoa genome assembly employing genomic variation for guided scaffolding

We propose to use the natural variation between individuals of a population for genome assembly scaffolding. In today's genome projects, multiple accessions get sequenced, leading to variant catalogs. Using such information to improve genome assemblies is attractive both cost-wise as well as scientifically, because the value of an assembly increases with its contiguity. We conclude that haplotype information is a valuable resource to group and order contigs toward the generation of pseudomolecules. Quinoa (Chenopodium quinoa) has been under cultivation in Latin America for more than 7500 years. Recently, quinoa has gained increasing attention due to its stress resistance and its nutritional value. We generated a novel quinoa genome assembly for the Bolivian accession CHEN125 using PacBio long-read sequencing data (assembly size 1.32 Gbp, initial N50 size 608 kbp). Next, we re-sequenced 50 quinoa accessions from Peru and Bolivia. This set of accessions differed at 4.4 million single-nucleotide variant (SNV) positions compared to CHEN125 (1.4 million SNV positions on average per accession). We show how to exploit variation in accessions that are distantly related to establish a genome-wide ordered set of contigs for guided scaffolding of a reference assembly. The method is based on detecting shared haplotypes and their expected continuity throughout the genome (i.e., the effect of linkage disequilibrium), as an extension of what is expected in mapping populations where only a few haplotypes are present. We test the approach using Arabidopsis thaliana data from different populations. After applying the method on our CHEN125 quinoa assembly we validated the results with mate-pairs, genetic markers, and another quinoa assembly originating from a Chilean cultivar. We show consistency between these information sources and the haplotype-based relations as determined by us and obtain an improved assembly with an N50 size of 1079 kbp and ordered contig groups of up to 39.7 Mbp. We conclude that haplotype information in distantly related individuals of the same species is a valuable resource to group and order contigs according to their adjacency in the genome toward the generation of pseudomolecules.

PhylomeDB V5: an expanding repository for genome-wide catalogues of annotated gene phylogenies

PhylomeDB is a unique knowledge base providing public access to minable and browsable catalogues of pre-computed genome-wide collections of annotated sequences, alignments and phylogenies (i.e. phylomes) of homologous genes, as well as to their corresponding phylogeny-based orthology and paralogy relationships. In addition, PhylomeDB trees and alignments can be downloaded for further processing to detect and date gene duplication events, infer past events of inter-species hybridization and horizontal gene transfer, as well as to uncover footprints of selection, introgression, gene conversion, or other relevant evolutionary processes in the genes and organisms of interest. Here, we describe the latest evolution of PhylomeDB (version 5). This new version includes a newly implemented web interface and several new functionalities such as optimized searching procedures, the possibility to create user-defined phylome collections, and a fully redesigned data structure. This release also represents a significant core data expansion, with the database providing access to 534 phylomes, comprising over 8 million trees, and homology relationships for genes in over 6000 species. This makes PhylomeDB the largest and most comprehensive public repository of gene phylogenies. PhylomeDB is available at http://www.phylomedb.org.

FA-nf: A Functional Annotation Pipeline for Proteins from Non-Model Organisms Implemented in Nextflow

Functional annotation allows adding biologically relevant information to predicted features in genomic sequences, and it is, therefore, an important procedure of any de novo genome sequencing project. It is also useful for proofreading and improving gene structural annotation. Here, we introduce FA-nf, a pipeline implemented in Nextflow, a versatile computational workflow management engine. The pipeline integrates different annotation approaches, such as NCBI BLAST+, DIAMOND, InterProScan, and KEGG. It starts from a protein sequence FASTA file and, optionally, a structural annotation file in GFF format, and produces several files, such as GO assignments, output summaries of the abovementioned programs and final annotation reports. The pipeline can be broken easily into smaller processes for the purpose of parallelization and easily deployed in a Linux computational environment, thanks to software containerization, thus helping to ensure full reproducibility.

Identification of ABC transporter G subfamily in white lupin and functional characterization of L.albABGC29 in phosphorus use

BACKGROUND

White lupin (Lupinus albus) is a leguminous crop with elite adaptive ability in phosphorus-deficient soil and used as a model plant for studying phosphorus (P) use. However, the genetic basis of its adaptation to low P (LP) remains unclear. ATPase binding cassette (ABC) transports G subfamily play a crucial role in the transportation of biological molecules across the membrane. To date, identification of this subfamily has been analyzed in some plants, but no systematic analysis of these transporters in phosphorus acquisition is available for white lupin.

RESULTS

This study identified 66 ABCG gene family members in the white lupin genome using comprehensive approaches. Phylogenetic analysis of white lupin ABCG transporters revealed six subclades based on their counterparts in Arabidopsis, displaying distinct gene structure and motif distribution in each cluster. Influences of the whole genome duplication on the evolution of L.albABCGs were investigated in detail. Segmental duplications appear to be the major driving force for the expansion of ABCGs in white lupin. Analysis of the Ka/Ks ratios indicated that the paralogs of the L.albABCG subfamily members principally underwent purifying selection. However, it was found that L.albABCG29 was a result of both tandem and segmental duplications. Overexpression of L.albABCG29 in white lupin hairy root enhanced P accumulation in cluster root under LP and improved plant growth. Histochemical GUS staining indicated that L.albABCG29 expression increased under LP in white lupin roots. Further, overexpression of L.albABCG29 in rice significantly improved P use under combined soil drying and LP by improving root growth associated with increased rhizosheath formation.

CONCLUSION

Through systematic and comprehensive genome-wide bioinformatics analysis, including conserved domain, gene structures, chromosomal distribution, phylogenetic relationships, and gene duplication analysis, the L.albABCG subfamily was identified in white lupin, and L.albABCG29 characterized in detail. In summary, our results provide deep insight into the characterization of the L.albABCG subfamily and the role of L.albABCG29 in improving P use.

Genome-wide in silico analysis of long intergenic non-coding RNAs from rice peduncles at the heading stage

Long intergenic non-coding RNAs (lincRNAs) belong to the category of long non-coding RNAs (lncRNAs), originated from intergenic regions, which do not code for proteins. LincRNAs perform prominent role in regulation of gene expression during plant development and stress response by directly interacting with DNA, RNA, or proteins, or triggering production of small RNA regulatory molecules. Here, we identified 2973 lincRNAs and investigated their expression dynamics during peduncle elongation in two Indian rice cultivars, Pokkali and Swarna, at the time of heading. Differential expression analysis revealed common and cultivar-specific expression patterns, which we utilized to infer the lincRNA candidates with potential involvement in peduncle elongation and panicle exsertion. Their putative targets were identified using in silico prediction methods followed by pathway mapping and literature-survey based functional analysis. Further, to infer the mechanism of action, we identified the lincRNAs which potentially act as miRNA precursors or target mimics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01059-2.

Genomic basis of high-altitude adaptation in Tibetan Prunus fruit trees

The Great Himalayan Mountains and their foothills are believed to be the place of origin and development of many plant species. The genetic basis of adaptation to high plateaus is a fascinating topic that is poorly understood at the population level. We comprehensively collected and sequenced 377 accessions of Prunus germplasm along altitude gradients ranging from 2,067 to 4,492 m in the Himalayas. We de novo assembled three high-quality genomes of Tibetan Prunus species. A comparative analysis of Prunus genomes indicated a remarkable expansion of the SINE retrotransposons occurred in the genomes of Tibetan species. We observed genetic differentiation between Tibetan peaches from high and low altitudes and that genes associated with light stress signaling, especially UV stress signaling, were enriched in the differentiated regions. By profiling the metabolomes of Tibetan peach fruit, we determined 379 metabolites had significant genetic correlations with altitudes and that in particular phenylpropanoids were positively correlated with altitudes. We identified 62 Tibetan peach-specific SINEs that colocalized with metabolites differentially accumualted in Tibetan relative to cultivated peach. We demonstrated that two SINEs were inserted in a locus controlling the accumulation of 3-O-feruloyl quinic acid. SINE1 was specific to Tibetan peach. SINE2 was predominant in high altitudes and associated with the accumulation of 3-O-feruloyl quinic acid. These genomic and metabolic data for Prunus populations native to the Himalayan region indicate that the expansion of SINE retrotransposons helped Tibetan Prunus species adapt to the harsh environment of the Himalayan plateau by promoting the accumulation of beneficial metabolites.

Identification and Expression Analysis of the Genes Involved in the Raffinose Family Oligosaccharides Pathway of Phaseolus vulgaris and Glycine max

Raffinose family oligosaccharides (RFO) play an important role in plants but are also considered to be antinutritional factors. A profound understanding of the galactinol and RFO biosynthetic gene families and the expression patterns of the individual genes is a prerequisite for the sustainable reduction of the RFO content in the seeds, without compromising normal plant development and functioning. In this paper, an overview of the annotation and genetic structure of all galactinol- and RFO biosynthesis genes is given for soybean and common bean. In common bean, three galactinol synthase genes, two raffinose synthase genes and one stachyose synthase gene were identified for the first time. To discover the expression patterns of these genes in different tissues, two expression atlases have been created through re-analysis of publicly available RNA-seq data. De novo expression analysis through an RNA-seq study during seed development of three varieties of common bean gave more insight into the expression patterns of these genes during the seed development. The results of the expression analysis suggest that different classes of galactinol- and RFO synthase genes have tissue-specific expression patterns in soybean and common bean. With the obtained knowledge, important galactinol- and RFO synthase genes that specifically play a key role in the accumulation of RFOs in the seeds are identified. These candidate genes may play a pivotal role in reducing the RFO content in the seeds of important legumes which could improve the nutritional quality of these beans and would solve the discomforts associated with their consumption.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

Control of the Rhizobia Nitrogen-Fixing Symbiosis by Common Bean MADS-Domain/AGL Transcription Factors

Plants MADS-domain/AGL proteins constitute a large transcription factor (TF) family that controls the development of almost every plant organ. We performed a phylogeny of (ca. 500) MADS-domain proteins from Arabidopsis and four legume species. We identified clades with Arabidopsis MADS-domain proteins known to participate in root development that grouped legume MADS-proteins with similar high expression in roots and nodules. In this work, we analyzed the role of AGL transcription factors in the common bean (Phaseolus vulgaris) - Rhizobium etli N-fixing symbiosis. Sixteen P. vulgaris AGL genes (PvAGL), out of 93 family members, are expressed - at different levels - in roots and nodules. From there, we selected the PvAGL gene denominated PvFUL-like for overexpression or silencing in composite plants, with transgenic roots and nodules, that were used for phenotypic analysis upon inoculation with Rhizobium etli. Because of sequence identity in the DNA sequence used for RNAi-FUL-like construct, roots, and nodules expressing this construct -referred to as RNAi_AGL- showed lower expression of other five PvAGL genes highly expressed in roots/nodules. Contrasting with PvFUL-like overexpressing plants, rhizobia-inoculated plants expressing the RNAi_AGL silencing construct presented affection in the generation and growth of transgenic roots from composite plants, both under non-inoculated or rhizobia-inoculated condition. Furthermore, the rhizobia-inoculated plants showed decreased rhizobial infection concomitant with the lower expression level of early symbiotic genes and increased number of small, ineffective nodules that indicate an alteration in the autoregulation of the nodulation symbiotic process. We propose that the positive effects of PvAGL TF in the rhizobia symbiotic processes result from its potential interplay with NIN, the master symbiotic TF regulator, that showed a CArG-box consensus DNA sequence recognized for DNA binding of AGL TF and presented an increased or decreased expression level in roots from non-inoculated plants transformed with OE_FUL or RNAi_AGL construct, respectively. Our work contributes to defining novel transcriptional regulators for the common bean - rhizobia N-fixing symbiosis, a relevant process for sustainable agriculture.

Genome-Wide Association Study and Genomic Prediction for Soybean Cyst Nematode Resistance in USDA Common Bean (Phaseolus vulgaris) Core Collection

Soybean cyst nematode (SCN, Heterodera glycines) has become the major yield-limiting biological factor in soybean production. Common bean is also a good host of SCN, and its production is challenged by this emerging pest in many regions such as the upper Midwest USA. The use of host genetic resistance has been the most effective and environmentally friendly method to manage SCN. The objectives of this study were to evaluate the SCN resistance in the USDA common bean core collection and conduct a genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers with SCN resistance. A total of 315 accessions of the USDA common bean core collection were evaluated for resistance to SCN HG Type 0 (race 6). The common bean core set was genotyped with the BARCBean6K_3 Infinium BeadChips, consisting of 4,654 SNPs. Results showed that 15 accessions were resistant to SCN with a Female Index (FI) at 4.8 to 9.4, and 62 accessions were moderately resistant (10 < FI < 30) to HG Type 0. The association study showed that 11 SNP markers, located on chromosomes Pv04, 07, 09, and 11, were strongly associated with resistance to HG Type 0. GWAS was also conducted for resistance to HG Type 2.5.7 and HG Type 1.2.3.5.6.7 based on the public dataset (N = 276), consisting of a diverse set of common bean accessions genotyped with the BARCBean6K_3 chip. Six SNPs associated with HG Type 2.5.7 resistance on Pv 01, 02, 03, and 07, and 12 SNPs with HG Type 1.2.3.5.6.7 resistance on Pv 01, 03, 06, 07, 09, 10, and 11 were detected. The accuracy of genomic prediction (GP) was 0.36 to 0.49 for resistance to the three SCN HG types, indicating that genomic selection (GS) of SCN resistance is feasible. This study provides basic information for developing SCN-resistant common bean cultivars, using the USDA core germ plasm accessions. The SNP markers can be used in molecular breeding in common beans through marker-assisted selection (MAS) and GS.

Gene Flow in Phaseolus Beans and Its Role as a Plausible Driver of Ecological Fitness and Expansion of Cultigens

The genus Phaseolus, native to the Americas, is composed of more than eighty wild species, five of which were domesticated in pre-Columbian times. Since the beginning of domestication events in this genus, ample opportunities for gene flow with wild relatives have existed. The present work reviews the extent of gene flow in the genus Phaseolus in primary and secondary areas of domestication with the aim of illustrating how this evolutionary force may have conditioned ecological fitness and the widespread adoption of cultigens. We focus on the biological bases of gene flow in the genus Phaseolus from a spatial and time perspective, the dynamics of wild-weedy-crop complexes in the common bean and the Lima bean, the two most important domesticated species of the genus, and the usefulness of genomic tools to detect inter and intraspecific introgression events. In this review we discuss the reproductive strategies of several Phaseolus species, the factors that may favor outcrossing rates and evidence suggesting that interspecific gene flow may increase ecological fitness of wild populations. We also show that wild-weedy-crop complexes generate genetic diversity over which farmers are able to select and expand their cultigens outside primary areas of domestication. Ultimately, we argue that more studies are needed on the reproductive biology of the genus Phaseolus since for most species breeding systems are largely unknown. We also argue that there is an urgent need to preserve wild-weedy-crop complexes and characterize the genetic diversity generated by them, in particular the genome-wide effects of introgressions and their value for breeding programs. Recent technological advances in genomics, coupled with agronomic characterizations, may make a large contribution.

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

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