Natural variation and gene regulatory basis for the responses of asparagus beans to soil drought

Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

Crop Landraces and Indigenous Varieties: A Valuable Source of Genes for Plant Breeding

Landraces and indigenous varieties comprise valuable sources of crop species diversity. Their utilization in plant breeding may lead to increased yield and enhanced quality traits, as well as resilience to various abiotic and biotic stresses. Recently, new approaches based on the rapid advancement of genomic technologies such as deciphering of pangenomes, multi-omics tools, marker-assisted selection (MAS), genome-wide association studies (GWAS), and CRISPR/Cas9 gene editing greatly facilitated the exploitation of landraces in modern plant breeding. In this paper, we present a comprehensive overview of the implementation of new genomic technologies and highlight their importance in pinpointing the genetic basis of desirable traits in landraces and indigenous varieties of annual, perennial herbaceous, and woody crop species cultivated in the Mediterranean region. The need for further employment of advanced -omic technologies to unravel the full potential of landraces and indigenous varieties underutilized genetic diversity is also indicated. Ultimately, the large amount of genomic data emerging from the investigation of landraces and indigenous varieties reveals their potential as a source of valuable genes and traits for breeding. The role of landraces and indigenous varieties in mitigating the ongoing risks posed by climate change in agriculture and food security is also highlighted.

Genome-wide association analysis reveals the optimal genomic regions for pod size in bean

The snap bean is the most commonly grown vegetable legume worldwide, and its pod size is both an important yield and appearance quality trait. However, the improvement of pod size in snap beans grown in China has been largely hindered by a lack of information on the specific genes that determine pod size. In this study, we identified 88 snap bean accessions and evaluated their pod size traits. Through a genome-wide association study (GWAS), 57 single nucleotide polymorphisms (SNPs) significantly associated with pod size were detected. Candidate gene analysis showed that cytochrome P450 family genes, WRKY, and MYB transcription factors were the predominant candidate genes for pod development, and eight of these 26 candidate genes showed relatively higher expression patterns in flowers and young pods. A significant pod length (PL) SNP and a single pod weight (SPW) SNP were successfully converted into kompetitive allele-specific polymerase chain reaction (KASP) markers and validated in the panel. These results enhance our understanding of the genetic basis of pod size, and also provide genetic resources for the molecular breeding of pod size in snap beans.

Understanding water conservation vs. profligation traits in vegetable legumes through a physio-transcriptomic-functional approach

Vegetable soybean and cowpea are related warm-season legumes showing contrasting leaf water use behaviors under similar root drought stresses, whose mechanisms are not well understood. Here we conducted an integrative phenomic-transcriptomic study on the two crops grown in a feedback irrigation system that enabled precise control of soil water contents. Continuous transpiration rate monitoring demonstrated that cowpea used water more conservatively under earlier soil drought stages, but tended to maintain higher transpiration under prolonged drought. Interestingly, we observed a soybean-specific transpiration rate increase accompanied by phase shift under moderate soil drought. Time-series transcriptomic analysis suggested a dehydration avoidance mechanism of cowpea at early soil drought stage, in which the VuHAI3 and VuTIP2;3 genes were suggested to be involved. Multifactorial gene clustering analysis revealed different responsiveness of genes to drought, time of day and their interactions between the two crops, which involved species-dependent regulation of the circadian clock genes. Gene network analysis identified two co-expression modules each associated with transpiration rate in cowpea and soybean, including a pair of negatively correlated modules between species. Module hub genes, including the ABA-degrading gene GmCYP707A4 and the trehalose-phosphatase/synthase gene VuTPS9 were identified. Inter-modular network analysis revealed putative co-players of the hub genes. Transgenic analyses verified the role of VuTPS9 in regulating transpiration rate under osmotic stresses. These findings propose that species-specific transcriptomic reprograming in leaves of the two crops suffering similar soil drought was not only a result of the different drought resistance level, but a cause of it.

Genome and pan-genome assembly of asparagus bean (Vigna unguiculata ssp. sesquipedialis) reveal the genetic basis of cold adaptation

Asparagus bean (Vigna unguiculata ssp. sesquipedialis) is an important cowpea subspecies. We assembled the genomes of Ningjiang 3 (NJ, 550.31 Mb) and Dubai bean (DB, 564.12 Mb) for comparative genomics analysis. The whole-genome duplication events of DB and NJ occurred at 64.55 and 64.81 Mya, respectively, while the divergence between soybean and Vigna occurred in the Paleogene period. NJ genes underwent positive selection and amplification in response to temperature and abiotic stress. In species-specific gene families, NJ is mainly enriched in response to abiotic stress, while DB is primarily enriched in respiration and photosynthesis. We established the pan-genomes of four accessions (NJ, DB, IT97K-499-35 and Xiabao II) and identified 20,336 (70.5%) core genes present in all the accessions, 6,507 (55.56%) variable genes in two individuals, and 2,004 (6.95%) unique genes. The final pan genome is 616.35 Mb, and the core genome is 399.78 Mb. The variable genes are manifested mainly in stress response functions, ABC transporters, seed storage, and dormancy control. In the pan-genome sequence variation analysis, genes affected by presence/absence variants were enriched in biological processes associated with defense responses, immune system processes, signal transduction, and agronomic traits. The results of the present study provide genetic data that could facilitate efficient asparagus bean genetic improvement, especially in producing cold-adapted asparagus bean.

Functional physiological phenotyping and transcriptome analysis provide new insight into strawberry growth and water consumption

Global warming is expected to increase agricultural water scarcity; thus, optimized irrigation schedules are important and timely for sustainable crop production. Deficit irrigation, which balances crop growth and water consumption, has been proposed, but the critical threshold is not easily quantified. Here, we conducted experiments on strawberry plants subjecting progressive drought following various water recovery treatments on the high-throughput physiological phenotyping system "Plantarray". The critical soil water contents (θcri), below which the plant transpiration significantly decreased, were calculated from the inflection point of the transpiration rate (Tr) - volumetric soil water content (VWC) curve fitted by a piecewise function. The physiological traits of water relations were compared between the well-watered plants (CK), plants subjecting the treatment of rewatering at the point of θcri following progressive drought (WR_θcri), and the plants subjecting the treatment of rewatering at severe drought following progressive drought (WR_SD). The results showed that midday Tr, daily transpiration (E), and biomass gain of the plants under WR_θcri treatment were equivalent to CK during the whole course of the experiment, but those under WR_SD treatment were significantly lower than CK during the water stress phase that could not recover even after rehydration. To explore the gene regulatory mechanisms, transcriptome analysis of the samples collected 12 h before, 12 h post and 36 h post water recovery in the three treatments was conducted. GO and KEGG enrichment analyses for the differentially expressed genes indicated that genes involved in mineral absorption and flavonoid biosynthesis were among the most striking transcriptionally reversible genes under the WR_θcri treatment. Functional physiological phenotyping and transcriptome data provide new insight into a potential, quantitative, and balanceable water-saving strategy for strawberry irrigation and other agricultural crops.

Unraveling the Genetic Architecture of Two Complex, Stomata-Related Drought-Responsive Traits by High-Throughput Physiological Phenotyping and GWAS in Cowpea (Vigna. Unguiculata L. Walp)

Drought is one of the most devasting and frequent abiotic stresses in agriculture. While many morphological, biochemical and physiological indicators are being used to quantify plant drought responses, stomatal control, and hence the transpiration and photosynthesis regulation through it, is of particular importance in marking the plant capacity of balancing stress response and yield. Due to the difficulties in simultaneous, large-scale measurement of stomatal traits such as sensitivity and speed of stomatal closure under progressive soil drought, forward genetic mapping of these important behaviors has long been unavailable. The recent emerging phenomic technologies offer solutions to identify the water relations of whole plant and assay the stomatal regulation in a dynamic process at the population level. Here, we report high-throughput physiological phenotyping of water relations of 106 cowpea accessions under progressive drought stress, which, in combination of genome-wide association study (GWAS), enables genetic mapping of the complex, stomata-related drought responsive traits “critical soil water content” (θ_cri) and “slope of transpiration rate declining” (K_Tr). The 106 accessions showed large variations in θ_cri and K_Tr, indicating that they had broad spectrum of stomatal control in response to soil water deficit, which may confer them different levels of drought tolerance. Univariate GWAS identified six and fourteen significant SNPs associated with θ_cri and K_Tr, respectively. The detected SNPs distributed in nine chromosomes and accounted for 8.7–21% of the phenotypic variation, suggesting that both stomatal sensitivity to soil drought and the speed of stomatal closure to completion were controlled by multiple genes with moderate effects. Multivariate GWAS detected ten more significant SNPs in addition to confirming eight of the twenty SNPs as detected by univariate GWAS. Integrated, a final set of 30 significant SNPs associated with stomatal closure were reported. Taken together, our work, by combining phenomics and genetics, enables forward genetic mapping of the genetic architecture of stomatal traits related to drought tolerance, which not only provides a basis for molecular breeding of drought resistant cultivars of cowpea, but offers a new methodology to explore the genetic determinants of water budgeting in crops under stressful conditions in the phenomics era.

Transcriptomic analysis of salt tolerance-associated genes and diversity analysis using indel markers in yardlong bean (Vigna unguiculata ssp. sesquipedialis)

BACKGROUND

High salinity is a devastating abiotic stresses for crops. To understand the molecular basis of salinity stress in yardlong bean (Vigna unguiculata ssp. sesquipedalis), and to develop robust markers for improving this trait in germplasm, whole transcriptome RNA sequencing (RNA-seq) was conducted to compare the salt-tolerant variety Suzi 41 and salt-sensitive variety Sujiang 1419 under normal and salt stress conditions.

RESULTS

Compared with controls, 417 differentially expressed genes (DEGs) were identified under exposure to high salinity, including 42 up- and 11 down-regulated DEGs in salt-tolerant Suzi 41 and 186 up- and 197 down-regulated genes in salt-sensitive Sujiang 1419, validated by qRT-PCR. DEGs were enriched in "Glycolysis/Gluconeogenesis" (ko00010), "Cutin, suberine and wax biosynthesis" (ko00073), and "phenylpropanoid biosynthesis" (ko00940) in Sujiang 1419, although "cysteine/methionine metabolism" (ko00270) was the only pathway significantly enriched in salt-tolerant Suzi 41. Notably, AP2/ERF, LR48, WRKY, and bHLH family transcription factors (TFs) were up-regulated under high salt conditions. Genetic diversity analysis of 84 yardlong bean accessions using 26 InDel markers developed here could distinguish salt-tolerant and salt-sensitive varieties.

CONCLUSIONS

These findings show a limited set of DEGs, primarily TFs, respond to salinity stress in V. unguiculata, and that these InDels associated with salt-inducible loci are reliable for diversity analysis.

Functional physiological phenotyping with functional mapping: A general framework to bridge the phenotype-genotype gap in plant physiology

The recent years have witnessed the emergence of high-throughput phenotyping techniques. In particular, these techniques can characterize a comprehensive landscape of physiological traits of plants responding to dynamic changes in the environment. These innovations, along with the next-generation genomic technologies, have brought plant science into the big-data era. However, a general framework that links multifaceted physiological traits to DNA variants is still lacking. Here, we developed a general framework that integrates functional physiological phenotyping (FPP) with functional mapping (FM). This integration, implemented with high-dimensional statistical reasoning, can aid in our understanding of how genotype is translated toward phenotype. As a demonstration of method, we implemented the transpiration and soil-plant-atmosphere measurements of a tomato introgression line population into the FPP-FM framework, facilitating the identification of quantitative trait loci (QTLs) that mediate the spatiotemporal change of transpiration rate and the test of how these QTLs control, through their interaction networks, phenotypic plasticity under drought stress.

Developments and Prospects in Imperative Underexploited Vegetable Legumes Breeding: A Review

Vegetable legumes are an essential source of carbohydrates, vitamins, and minerals, along with health-promoting bioactive chemicals. The demand for the use of either fresh or processed vegetable legumes is continually expanding on account of the growing consumer awareness about their well-balanced diet. Therefore, sustaining optimum yields of vegetable legumes is extremely important. Here we seek to present d etails of prospects of underexploited vegetable legumes for food availability, accessibility, and improved livelihood utilization. So far research attention was mainly focused on pulse legumes' performance as compared to vegetable legumes. Wild and cultivated vegetable legumes vary morphologically across diverse habitats. This could make them less known, underutilized, and underexploited, and make them a promising potential nutritional source in developing nations where malnutrition still exists. Research efforts are required to promote underexploited vegetable legumes, for improving their use to feed the ever-increasing population in the future. In view of all the above points, here we have discussed underexploited vegetable legumes with tremendous potential; namely, vegetable pigeon pea (Cajanus cajan), cluster bean (Cyamopsis tetragonoloba), winged bean (Psophocarpus tetragonolobus), dolichos bean (Lablab purpureus), and cowpea (Vigna unguiculata), thereby covering the progress related to various aspects such as pre-breeding, molecular markers, quantitative trait locus (QTLs), genomics, and genetic engineering. Overall, this review has summarized the information related to advancements in the breeding of vegetable legumes which will ultimately help in ensuring food and nutritional security in developing nations.

Orphan genes are involved in drought adaptations and ecoclimatic-oriented selections in domesticated cowpea

Orphan genes (OGs) are genes that are restricted to a single species or a particular taxonomic group. To date, little is known about the functions of OGs in domesticated crops. Here, we report our findings on the relationships between OGs and environmental adaptation in cowpea (Vigna unguiculata). We identified 578 expressed OGs, of which 73.2% were predicted to be non-coding. Transcriptomic analyses revealed a high rate of OGs that were drought inducible in roots when compared with conserved genes. Co-expression analysis further revealed the possible involvement of OGs in stress response pathways. Overexpression of UP12_8740, a drought-inducible OG, conferred enhanced tolerance to osmotic stresses and soil drought. By combining Capture-Seq and fluorescence-based Kompetitive allele-specific PCR (KASP), we efficiently genotyped single nucleotide polymorphisms (SNPs) on OGs across a 223 accession cowpea germplasm collection. Population genomic parameters, including polymorphism information content (PIC), expected heterozygosity (He), nucleotide diversity (π), and Tajima's D statistics, that were calculated based on these SNPs, showed distinct signatures between the grain- and vegetable-type subpopulations of cowpea. This study reinforces the idea that OGs are a valuable resource for identifying new genes related to species-specific environmental adaptations and fosters new insights that artificial selection on OGs might have contributed to balancing the adaptive and agronomic traits in domesticated crops in various ecoclimatic conditions.

Variation in relative water content, proline accumulation and stress gene expression in two cowpea landraces under drought

Many landraces of cowpea [Vigna unguiculata (L.) Walp.] are adapted to particular geographical and climatic conditions. Here we describe two landraces grown respectively in arid and temperate areas of Algeria and assess their physiological and molecular responses to drought stress. As expected, when deprived of water cowpea plants lose water over time with a gradual reduction in transpiration rate. The landraces differed in their relative water content (RWC) and whole plant transpiration rate. The landrace from Menia, an arid area, retained more water in adult leaves. Both landraces responded to drought stress at the molecular level by increasing expression of stress-related genes in aerial parts, including proline metabolism genes. Expression of gene(s) encoding proline synthesis enzyme P5CS was up regulated and gene expression of ProDH, a proline catabolism enzyme, was down regulated. Relatively low amounts of proline accumulated in adult leaves with slight differences between the two landraces. During drought stress the most apical part of plants stayed relatively turgid with a high RWC compared to distal parts that wilted. Expression of key stress genes was higher and more proline accumulated at the apex than in distal leaves indicating that cowpea has a non-uniform stress response at the whole plant level. Our study reveals a developmental control of water stress through preferential proline accumulation in the upper tier of the cowpea plant. We also conclude that cowpea landraces display physiological adaptations to water stress suited to the arid and temperate climates in which they are cultivated.

Genomic regions, cellular components and gene regulatory basis underlying pod length variations in cowpea (V. unguiculata L. Walp)

Cowpea (V. unguiculata L. Walp) is a climate resilient legume crop important for food security. Cultivated cowpea (V. unguiculata L) generally comprises the bushy, short-podded grain cowpea dominant in Africa and the climbing, long-podded vegetable cowpea popular in Asia. How selection has contributed to the diversification of the two types of cowpea remains largely unknown. In the current study, a novel genotyping assay for over 50 000 SNPs was employed to delineate genomic regions governing pod length. Major, minor and epistatic QTLs were identified through QTL mapping. Seventy-two SNPs associated with pod length were detected by genome-wide association studies (GWAS). Population stratification analysis revealed subdivision among a cowpea germplasm collection consisting of 299 accessions, which is consistent with pod length groups. Genomic scan for selective signals suggested that domestication of vegetable cowpea was accompanied by selection of multiple traits including pod length, while the further improvement process was featured by selection of pod length primarily. Pod growth kinetics assay demonstrated that more durable cell proliferation rather than cell elongation or enlargement was the main reason for longer pods. Transcriptomic analysis suggested the involvement of sugar, gibberellin and nutritional signalling in regulation of pod length. This study establishes the basis for map-based cloning of pod length genes in cowpea and for marker-assisted selection of this trait in breeding programmes.